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CN ' Abstract 

In this paper we give examples of applications of general methods 
, of quantization by symmetrization of classical integrable systems, which 

have been illustrated in two previous works by the same authors. We 
consider two classes of systems in n spatial dimensions, which respectively 
describe a point particle in a central force field and a freely rotating rigid 
body. In the former case, the application of the general methods to an 
integrable classical system leads in an almost straightforward way to the 
quasi-integrability of the corresponding quantum system. In the latter 
case instead, a modification of the symmetrization procedure is necessary 
in order to achieve quantum integrability for n = 6. 
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00 '■ 1 Introduction 

In two previous papers of this series, we have introduced the concept of quasi- 
• integrable quantum system [J, and we have established general methods to 

I obtain examples of such systems starting from classical integrable systems [5]. 

Integrability of an operator H is defined as the existence of a sufficiently large 
set F of operators which commute with H, more exactly a quasi-integrable set F 
. ■ of operators. The main source of integrable sets are Lie closed sets of operators 

r> I commuting with H. Making the union of several such sets one can obtain an 

. integrable set F. An important particular case of Lie closed set is a Lie algebra 

of operators. General methods for the constructions of integrable sets are based 
on the symmetrization of the products of operators, which correspond to the 
elements of an integrable set of functions for the classical system. In the present 
paper these methods will be applied to some important classes of integrable 
classical and quantum systems. 

In section[5]we consider systems in a Euclidean space of arbitrary dimension 
n, which describe a point particle in a central force field. The discussion is then 
extended to more general one-particle systems which are symmetric with respect 
to the group of rotations SO{n). For all these cases we construct various types 
of classical integrable sets of functions. These sets contain in general 2n — k 
elements, where k is equal to the number of elements in the central subset 
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[3]. This number, for the various integral sets here considered, can take aU 
possible values from 2 to n. Each integrable set can be applied to all systems 
whose hamiltonian is an arbitrary function of the central elements. In all these 
cases we show that, by applying the general results of one can obtain a 
corresponding integrable quantum systems with an equal number k of central 
operators. 

In section [3] we then consider a freely rotating rigid body in a Euclidean 
space of arbitrary dimension n. In the classical case, it is known that this 
system is completely integrable. By applying our scheme of noncommutative 
integrability [3l SJ [5] , we find how the number k of central integrals depends on 
the space dimension n and on the properties of the set of generalized moments 
of inertia of the body. We then show that, for n < 5, the application of the 
general results of [5] leads in an almost straightforward way to the integrability 
of the corresponding quantum systems. However, for n — 6 we find that a 
modification of the symmetrization procedure is necessary in order to obtain 
a quasi-integrable set of operators. In fact, a Manakov polynomial of fourth 
order in the left-invariant momenta, which belongs to the central subset of 
the classical integrable set, does not commute with the hamiltonian operator 
after symmetrization in the momenta. However, commutativity can be restored 
by adding to it a suitable second order polynomial. One can conjecture that 
analogous procedures can be applied also for n > 6. 

2 One-particle systems in a n-dimensional cen- 
tral force field 

It is well-known that several mechanical systems have been proven to be in- 
tegrable both at the classical and at the quantum level. The integrability of 
various classes of systems is discussed for instance in [HI [3 |H1 [SI [13 [HI [H]- 
These systems usually consist of point particles moving in a space of one or 
more dimensions, subjected to an external potential or mutually interacting via 
suitable two-particle potentials. In particular, in |13l a class of maximally super- 
integrable systems is studied, which includes as a particular case the hydrogen 
atom in n dimensions. In this section we too shall consider systems in n di- 
mensions which generalize in some sense the hydrogen atom problem, although 
our aim will be partly different with respect to most of the cited investigations. 
We shall not in fact restrict ourselves to considering hamiltonians which are the 
sum of the usual kinetic term and of a potential term dependent only on the 
position. We shall instead consider a generic invariant hamiltonian with respect 
to the group of n-dimensional rotations, and we shall look for all the possible 
integrable sets of functions and operators which can be constructed for such an 
hamiltonian. It is then obvious that each of these sets can also be associated 
with the entire class of integrable systems, whose hamiltonian is expressible as 
a functions of the central elements of the set. In this way our approach leads 
to the systematic individuation of families of integrable systems in n spatial 
dimensions. However, we shall not discuss the possible physical interpretation 
of the systems obtained with this method. 
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2.1 Classical particle in a central force field 

The hamiltonian function of this system in an n-dimensional euclidean space 
has the foUowing form: 

H=^p^ + U{r) (2.1) 

where U e C°°(0, +cx)). Here we use the notation 

r = \/^, X = (xi, . . . ,x„) , p = {pi, . . . ,pn) . 

The first term on the right-hand side of (I2.ip corresponds to the kinetic energy 
of the particle, and the second one corresponds to its potential energy. The con- 
figuration space K of this classical system is the n-dimensional euclidean linear 
space M", more exactly, K = R"\{0}. The group of orthogonal transformations 
G = SO{n) acts on this space. This action in an orthonormal basis in E" is de- 
fined by orthogonal matrices. The dimension of this Lie group is = n{n—l)/2. 
The action of this Lie group G transfers onto the cotangent bundle T*R" — 
to R". This bundle, without the cotangent space Tq*R" to R" at the point 0, 
represents the phase space M = T*K of the classical system. The action of this 
group G conserves the hamiltonian function H = H{x,p) = ^p^ + U{r) of the 
classical system. 

Let us denote the Lie algebra of the group G as g = so(n). Each element 
a of g is associated with a vector field Va on R". The corresponding vector 
field on the symplectic manifold M C T*R" is hamiltonian with hamiltonian 
function Pa (™) — {p,Va{x)), where p — TH ^ T^'M.'^ is a linear form on T^j^M^. Let 
{vai{x)^ . . . , Van{x)) bc the components of the vector field Va in coordinates x — 
); then Pa{p, x) = J^i PiVai{x), where (p, x) are canonical coordinates 
on r*R". In the considered case, in which the group is G = SO{n), an element 
a of the Lie algebra g is represented by a skew-symmetric matrix A = Aa , and 
the vector Va{x) at the point x = (xi, . . . , x„) has the form Va{x) — —AaX. The 
action of the group G on T*R" is a Poisson action, i.e., {Pa,Pb} = P[a,b]: where 
a, b are any two elements of g, and [a, b] is the commutator defined in this algebra. 
Each system of cartesian coordinates in R" defines in the algebra g a basis whose 
elements are given in these coordinates by matrices D^^ ^ ^ ^ i < j ^ n, having a 
particularly simple form. These matrices have only two non-zero elements which 
are equal to ±1. Namely, Dfj = —D^j- = 1, where D^-'; denotes the element of 
the matrix D^^ lying at the intersection of row k and column Making use of 
Kronecker symbol Sij, we can write in general 

-Dfe; = SkiSij - SkjSii . (2.2) 

According to this formula, a matrix D^^ can naturally be defined also for i > j, 
and we have D^^ ~ ~D^^ i,j — 1, . . . , n. The commutation relations between 
these matrices are 

^ Di^k^ = -5,hD^>' - J.feD*" + S^kD^'^ + SjhD''' . (2.3) 

Note that, whenever the commutator is nonzero, only one of the four terms on 
the right-hand side is different from zero. 

Let Pij denote the function Pa, corresponding to the matrix Aa = D^K Then 

Pij = XiPj - XjPi . (2.4) 
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From (|2.3p one immediately derives that the Poisson brackets relations between 
these functions are 

{Pij,Phk} = ~5ihPjk - SjkPih + SikPjh + SjhPik ■ (2.5) 

Let P{x,p) denote the iV-dimensional vector {Pij{x,p), I < i < j < n). It is 
easy to verify that the rank of the map P : — >■ at a typical point {x,p) 
is equal to 2n — 3. Let P^ denote the square of the length of vector P, i.e., 
:— J2i<i<j<n -^ij ~ ^^P^ ~ (-^ ' P)^- Here x ■ p denotes the scalar product of 
vectors x and p, i.e., x ■ p := X]r=i ^iPi- Using (|2.5p it is easy to check that 

{F2,P,,} = (2.6) 

for any component of vector P. It is clear that one can select 2n ~ A 
components L — {Pi^j^, ■ ■ ■ ,^12^-41271-4) of this vector, such that the set 11 :— 
{P^,L) defines a regular map 11 : — >• K^"~'^ almost everywhere in E^^. This 
means that the rank of the map 11 is equal to 2n — 3 almost everywhere. A 
possible choice is L = (P13, Pi4, ■ . ■ , Pi«, P23, -f24, • ■ . , P2n)- 

Let us add to this set the hamiltonian function H = H{x,p), and denote by 
F — {H, P^; L) the resulting set of 2n — 2 functions. It is easy to see that this 
set F is functionally independent almost everywhere, that is the set of critical 
points of the map defined by this set has zero measure, is a closed set and is 
nowhere dense. Using the relations 

{xi,Pjk} = SijXk - SikXj , (2.7) 
{Pt, Pjk} = SijPk - SikPj (2.8) 

for i, J, fc = 1, . . . , n, it is also easy to verify that 

{r^P}-0, {p^Pj^O, (2.9) 

that is {r^,Pjk} = {p^,Pjk} = 0. Hence {H,P} = 0. Since L d P, this implies 
that the set F has 2 central functions, H and P^. According to the definition 
given in |T], P is thus an integrable set with two central integrals, and the system 
with hamiltonian function H is globally integrable with set of invariants F. The 
conservation of P implies that the orbit of the particle lies in a 2-dimensional 
plane. 

Proposition 2.1. Let {Vi,...,Vi) and {Wi,...,Ws) be two sets of function- 
ally independent functions on a 2n- dimensional symplectic manifold, such that 
{Vi, Wk\ = for i — 1, . . . , / and k = 1, . . . , s. Then I + s < 2n. 

Taking into account the above well-known result, we can describe the set of 
all integrable classical systems which are invariant with respect to the action of 
the group G = SO{n) on M^p, i.e., the integrable systems whose hamiltonian H 
is in involution with vector P. 

Lemma 2.2. We have {H, P} — if and only if locally H — f{p'^,r, P^). 

Proof. Let us suppose that H — f{p'^,r,P^). Since the functions r,p'^,P^ are 
in involution with P (see above), we have {H, P} = 0. 

Viceversa, let us suppose that there exists a function H{x,p) such that 
{H,P} = and that has not locally the form H = f{p'^,r,P^). In this case 
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we would have 4 functionally independent functions p^, r, and H, which are 
in involution with the 2n — 3 functionally independent functions of the set 11. 
Since the sum of the numbers of functions belonging to these two sets equals 
271 + 1, this would be in contradiction with proposition □ 

Proposition 2.3. In the real analytic case, a hamiltonian function H is in 
involution with P and is integrable if and only if it has locally the following 
form: H ^ f{p^ ,r, P^). 

The situation considered in all the present article refers to the more gen- 
eral case of infinitely differentiable functions, that is of class C°° . In this 
case, if {H,P} — and H is integrable, then H — f{p^,r,P^). Viceversa, 
if H = f{p'^,r,P^), where the 2-covector df /dijP' ,r) — {df /d{p'^),df /dr) is 
not zero, more exactly df /d{p^,r) ^ almost everywhere, then H is integrable 
and {H,P} — 0. Such systems are always integrable with k = 2, with central 
integrals H and P^ . In certain cases it is possible to find an additional integral, 
and to have integrability with k = 1. For example in the case of the Newton 
potential, that is for H = p'^/2 — a/r, and in the case of identical uncoupled 
oscillators, that is for H = [p^ + r^)/2. The system with hamiltonian func- 
tion H which is a function only of the square angular momentum, more exactly 
H = f{P^), where df ^ Q almost everywhere, is integrable with k — 1. 

Let us present the integrable sets F = F(II) which correspond to these 
hamiltonian functions. For H = f(j)^,r,P'^), where the function f is not lo- 
cally functionally dependent only on P , one can take F = {H,P'^;L), and 
therefore k = 2. For H = p'^/2 — a/r one can take F — {H; P"^ , L, Ai), 
where Ai — 'YTi=i^^3'Pj ~ axi/r. For H = {p^ + r^)/2 one can take F — 
[H- Hi,H2,..., Hn-u Pi2, Pi3,---, Pin), where H, = [pj + a;|)/2. For H ^ 
f{P^) one can take F — {P^;p^,r,L). In all cases except the first one, we have 
k = l. 

Proof. According to the previous lemma, {H, P} = if and only if = 

If df /d{p^,r) ^ almost everywhere, we can repeat the proof given at the 
beginning of section \TJ\ for the case H = p'^ /2 — U{r). By adding the function 
H to the set 11 = (P^, L), we thus obtain a set of 2?t, — 2 functions which defines 
almost everywhere a regular map and has two central integrals, H and P^. 
Therefore the system is integrable with k = 2. If instead / — f{P'^), then the 
set F = r, L) is an integrable set with k — 1. 

For the Kepler system, H = p^/2 — a/r, it is straightforward to verify that 
{H,Ai} = 0, where 

n 

= ^ PijPj - a— = (p'^ - x^ - {x ■ p)pi 

for i — 1, . . . ,n J4j . Since vector A lies in the plane of the orbit and satisfies 
A^ = 2P^II+a^, obviously only one component of A is functionally independent 
of the set {H, P^, L). Hence F — (H; P^,L, Ai) is an integrable set with fc = 1. 
Of course, one has to keep in mind that the potential U{x) = —a/r is a solution 
of the Laplace equation Y^^=i d^U/dxj = only for n — i. 

The hamiltonian H — {p^ + r^)/2 actually describes a set of resonators with 
equal frequencies, and will be considered again in a following paper. 
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In all the considered cases, the linear independence of the differentials of 
the functions F{H) can be checked directly by considering the corresponding 
jacobian matrices. In the analytic case, either H — f{P^), or df /d{p^,r) ^ 
almost everywhere: therefore integrability is always guaranteed. □ 

Remark 2.1. In the previous proposition we have considered the potential a/r 
for a particle in a space of arbitrary dimension n. Of course one has to keep in 
mind that such a potential is a Green function for the n-dimensional Laplace 
operator only for n — 3. 

Let us consider the case H = f{p^,r,P^), with df/d{p'^,r) ^ almost 
everywhere. We have seen that these systems are integrable with k = 2. This 
means that the typical invariant surface for the phase-flow of the system is a 
two-dimensional torus. It is however possible to find for the same systems also 
integrable sets with a larger number k of central integrals, up to the maximum 
possible number k = n which corresponds to standard Liouville integrability. In 
this way one can construct a larger class of integrable systems, which includes 
all systems whose hamiltonian is an arbitrary function of the central elements 
of the set. In general, such systems will no longer be invariant under the action 
of the whole group SO{n), but only of some subgroup of it. 

For example, in the familiar case n — 3, one can take F = {H, P^; P12, P13), 
with fc = 2, but also F = {H, P^, F12), with fc 3, or in general F = {H, P'^,Pa), 
where Pa is the momentum associated with any arbitrary element a G so(3). 
It follows that any system with hamiltonian K — g{p^ ,r, P^ , Pa), where 5 is a 
function such that dg/d{p^,r) ^ 0, is integrable with fc = 3. Of course any 
such system is only invariant with respect to the one-parameter subgroup of the 
rotations around the axis associated with a. 

Definition 2.1. We say that two integrable sets are functionally equivalent if 
the elements of one set are locally functions of the elements of the other set. 

Owing to the arbitrariness in the choice of the element a € so(3), we see 
that, for the system with hamiltonian H = f{p^,r,P^), there exist infinitely 
many functionally inequivalent integrable sets with fc = 3. 

Let us now consider the case of arbitrary n. Note first of all that in the 
integrable set described in proposition 12. 3[ all n coordinates of configuration 
space are treated on the same footing, since all possible choices of the 2n — 4 
noncentral elements of the set actually lead to functionally equivalent integrable 
sets. More generally, under any transformation of the group SO{n), the set F 
with fc = 2 is transformed into an equivalent set. One can however construct 
other integrable sets in the following way. One takes the function P^ as central 
element, and then splits the set of n coordinates of configuration space into 
two arbitrary disjoint subsets. One takes as additional central elements the two 
functions, one for each of these two subsets, which are obtained by summing 
the squares of all the components of P acting on the coordinates of the subset. 
Then, for any of the two subsets of coordinates, one can proceed in two alter- 
native ways. Either one takes as integral functions a suitable set L' of 2n' — 4 
momenta acting on the coordinates of the subset, where n' is the number of such 
coordinates, or one splits again the subset into two arbitrary smaller disjoint 
subsets, and repeats the procedure. If one wishes, one can continue splitting the 
subsets into two parts, until one is left with only subsets consisting of either one 
or two space coordinates (the splitting of a set of two coordinates is ineffective 
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with respect to the resulting integrable set). Coming back to the case rt = 3, we 
see that the integrable system F = {H,P^,Pi2) considered above corresponds 
to the splitting of the set of coordinates {xi,X2,X3) into the two subsets {xi,X2) 
and {xs). 

The general procedure is described in a formal way by the following propo- 
sition. 

Proposition 2.4. For any n > 2 and for any z — 1, ... ,n — 1, it is possible 
to construct in a recursive manner sets Zn^z o,nd L„.z of polynomial functions 
of degree < 2 in the variables Pn '■= {Pij^l < i < j < n), with the following 
properties: 

1. Zn^z contains z elements, 

2. Ln^z contains 2{n — z — 1) elements, all of degree 1 in Pn, 

3. the set Yln.z '■— {Zn.z, Ln,z) is functionally independent, 

4. {z„,„n„,j = 0. 

Given any set A, it is useful to denote with (jA the number of its elements. 
Properties 1 and 2 can thus be written f,Zn^z = z and (|L„,z = 2{n — z — 1) 
respectively. Property 4 means that any element of Zn,z is in involution with 
all elements of the set Iln,z- 

Proof. According to proposition 12. 31 for z = 1 one can take Zn,i = Pn and form 
Ln.i by collecting 2n — 4 suitable elements of P„. For n = 2 we can only have 
z = 1, ^2,1 ~ {P12) and L2,i = 0. In order to construct sets Zn^z and Ln^z, with 
n > 2 and 2 < z < n — 1, we shall proceed by induction on n. 

Let us take m > 2 and suppose that, for all n = 2, . . . ,m — 1, we have 
constructed sets Zn,z and Ln.z of polynomial functions of degree < 2 in P„, for 
all possible z — — 1, satisfying properties 1-4 specified above. Let us 

split the set of indexes N„-i :— (1, . . . ,m) into two arbitrary nonempty disjoint 
subsets Ii and I2, such that jj/i = ni, jl/2 = n2, ni + n2 — m and Nm = 
Ii U l2- For k = 1,2, consider the two sets of momenta P^''^ C Pm, with 
PC^') := (Py , i,j e Ik, i < j) if 1 < nfc < m - 1, and P^'''> := if rife = 
1. This means that the elements of P*^'^' are the generators of the orthogonal 
transformations of the subspace R"'' C having set of coordinates {xi, i ^ Ik). 
We have obviously {p(^),p(^)} = 0. If 1 < < m — 1, consider for any 
Zfc, with 1 < Zfc < Hfe - 1, the sets Zn^^zAP'"^^) and Ln^^z^iP^'"''), which are 
obtained from the set of polynomials Zn^^zk and Ln^^zk by replacing the variables 
P„^ :— {Pij, I < i < j < nk) with P^''\ If instead Uk = 1, take Zk — 0, 
Zi^oiP^^')) and Li,o(P('''0 := 0. We thus have in all cases ttZ„,,2,(p('=)) = 
Zfc and tJi„fc,2fc(^'*''^) = 2(nfc - Zfc - 1). Finally, take 

Z = Zl + Z2 + 1 , 

z - (P^ Z CP^^h Z (P'-^^)) 
L - (L fP^^h L (P'^^^)) 

where P^ :— X]i<i<j<m Pfj- I* is easy to see that the sets Zm,z and Lm^z satisfy 
properties 1-4 above for n = ni. If rii = 1, zi = 0, 712 = m — 1, Z2 = 1, we 
obtain z = 2. With any other choice of rifc and z^, fc = 1, 2, z can assume any 
value from 3 to m — 1. □ 
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Let the hamiltonian of a system have the form H = f{p'^,r,P'^), with 
df I dijP' ,r) ^ almost everywhere. From proposition 12.41 it follows that the 
sets of functions Fn^z ■= {H, Zn^z', Ln,z) are integrable sets, with subset of cen- 
tral elements (iJ, Zn,z), for all z = 1, . . . , n — 1. We have tlFn,z — 2n — k and 
k — z + I. Hence the number k of central elements can take all values from 
k — 2 to k = n. 

It has to be noted that all integrable sets with k > 2, obtained by means 
of proposition I2.4[ depend on the choice of a cartesian set of coordinates on 
K". This means that, to any such set of coordinates, it corresponds in general 
an inequivalent integrable set. If one performs a transformation of SO{n) on 
configuration space, then a given integrable set is transformed into an equivalent 
one only if the transformation leaves invariant all the central functions of the 
set. 

As an example of application of proposition 12.41 in the two following tables 
we show explicitly some integrable sets which are obtained in the two cases 
71 = 4 and n = 5 respectively. In these tables we use the notation P('^]^23) 

_ p2 I p2 



P23 P(^i234) 



Ph 



P23 



P^i 



p2 



{H,P';Pi3,Pu,P23,P24) 

{H, P(2^23); -Pi2i -P13) 

{H, P'^ , P^^^^y P12) 



(iJ,p2,Pl2,P34) 

Table 1: Integrable sets for H = f{p^,r,P^) and n 



= 4. 



F 



{H, P^; Pi3, Pi4, ^15,^23, P24, P25) 

{H, P^, P^^^234); ^13, -P14, ^23, ^24) 
{H, P2 , P(^;^234) ' -^(^123) ; -^13 , ^23 ) 
{H,P ,-P(i234)''P'(123)''^12) 
{H, P2 , P(^i234) ' -^12 ' -^34) 
{H, P^ , P^i23) ' ^45 ; ^"12 , ^13 ) 
{H, P2 , P(^i23) ' ^45 , P12) 



Table 2: Integrable sets for H = f{p^,r,P^) and n = 5. 



2.2 Quantum particle in a central force field 

The hamiltonian operator H of this system is obtained from the hamiltonian 
function (|2.ip by standard quantization (see definition in [I]), which here simply 
consists in the substitution p — >• d/dx and the the replacement of the multiplica- 
tion of functions with the composition of corresponding operators, in symbols: 
X — > o. We thus obtain 

H=^f + U{r), (2.10) 



8 



where 



The operator -p^ is the Laplace operator in cartesian coordinates. We can pro- 
ceed as for the classical system, and on the basis of the classical formulas we will 
obtain the corresponding formulas where functions are converted into operators 
and Poisson brackets into Lie brackets. In a similar way we will also verify the 
quasi- independence of operators. Let us fix cartesian coordinates in M". Let 
(a;,p) be the corresponding canonical coordinates on = T*R". Let us con- 
sider the standard set of operators where p = d/dx. This set is obtained 
by canonical quantization from the coordinates (a;,p). Therefore, identifying 
(a;,p) and with the sets B and B respectively, one can study the algebra 
of the polynomial functions of these operators by applying propositions 3.1, 4.1, 
and also remark 4.4 of j2j. This remark is useful in order to deal with arbitrary 
functions of r. 

Let us consider the operators Py , 1 < i < j < n, obtained by symmetric 
quantization from the classical momenta Pij — XiPj — XjPi . Taking into account 
the canonical commutation relations 

[xi,Xj]=0, [pj,Pj]=0, [pi,Xj\^5ij (2.11) 

for i,j = 1, . . . , n, where 5ij is the Kronecker symbol, we have that = XiPj — 
XjPi, i.e., these operators coincide with the standard quantization of momenta 
Pij. From the quadratic dependence of Pij on {x,p), and from proposition 3.1 
of [2] (case 1), it follows that the commutation relations among the operators 
Pij have the same form as the Poisson brackets (|2.5p among the corresponding 
classical functions: 

[P^J , Phk] = ~6^hPjk - S,kP,h + 5,kPjh + SjhPk ■ (2.12) 

Similarly, from (|2.9p and proposition 4.1 of 2 it follows that 

[r2,P]=0, [p^P]=0. (2.13) 

It is also easy to verify that [U{r),P] = for any function U, in accordance with 
the first of ([^T^ and with remark 4.4 of [5]. We thus conclude that [H, P] = 0. 
Note that {p^, r^} = Ax ■ p ^ and correspondingly [p^,r'^] — 2{x ■ p + p ■ x) — 
4:X-p + 2nj^0. 

Using proposition 2.5 of [2 , it is easy to check that the operator (P^)''^™, 
obtained by symmetrization with respect to {x,p) of the square length P^ of 
momentum P, coincides with the operator P^ = J2i<j ^ij '^P additive 
constant. We have in fact P^ = (p2-jsym _|_ _ i)/4. Since the additive 
constant n{n — l)/4 is irrelevant for Lie brackets, from classical relations (|2.6p 
and from proposition 4.1 (case b) of [2^ we obtain 

[P^Py] = [(p2)^y-,Py] =0. (2.14) 

The quasi-independence of the set II = (P^,I/) of 2n — 3 operators follows 
from the functional independence of the corresponding set of symbols 11 = 
(P^, L), and from the homogeneity of these functions with respect to p. In fact. 
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the latter property implies that these functions coincide with their respective 
main parts. Let us consider the set F = {H,P^;L), where H is given by (|2.10p . 
If U{r) = 0, it is easy to check that the set of corresponding symbols F — 
{H,P^;L) is functionally independent. Since these functions are homogeneous 
with respect to p, these functions are quasi-independent. Furthermore, since 
the main part of H does not depend on U, the property of quasi-independence 
is true for arbitrary U. Therefore, the quantum system with hamiltonian H is 
quasi- integrable with integrable set F = {H, P^; L) and k = 2 central operators, 
H and 

It is possible to give a partial characterization of integrable quantum systems 
which are invariant with respect to the action of the group G = SO{n) on K^p, 
i.e., the systems whose hamiltonian operator H commutes with vector operator 
P. 

Proposition 2.5. If the hamiltonian operator H of a system has the form: H = 
f{p^, P^, (7i(r), . . . , gi(r)), where the function f is an arbitrary noncommutative 
polynomial in the 1 + 2 variables [p^ ,P'^ , gi{r), . . . , gi(r)), and gi(r), . . . , gi{r) 
are arbitrary functions ofr, then [H,P] = 0. 

Viceversa, let H be an arbitrary operator of class O on K = R" \ {0}, such 
that [H^P] = 0. Then the symbol of its main part MH with respect to linear 
momenta p (see definitions in fljl) has the form MH = g(ji'^,r,P^), where g 
is a homogeneous polynomial in the two variables {p^,P^), whose coefficients 
are arbitrary functions of r defined for all r > 0. If the polynomial g satisfies 
the condition dg/d{p^,r) ^ almost everywhere, then the systems is quasi- 
integrable with k — 2, with central integrals H and P^ . In certain cases it is 
possible to find an additional integral, and to have quasi-integrability with k = 1. 
This is the case for example for the Newton potential, that is for H = p^ /2 — a/r , 
and for identical uncoupled oscillators, that is for H = (p^ + r^)/2. The system 
with hamiltonian function H which is a function only of the square angular 
momentum, more exactly H = f{P'^), where df ^ almost everywhere, is 
quasi-integrable with k — I. 

Let us present the integrable sets of operators F — F{H) which corre- 
spond to integrable quantum systems with these hamiltonian operators. For 
H = f(p^,P^,gi{r),...,gi{r)), where the function f is not locally functionally 
dependent only on P^ , one can take F = {II,P^;L), and therefore k — 2. For 
H = p'^/2 — a/r one can take F = {H; P^ , L, Ai), where A\ — X]J=2(-^1jPj + 
PjPij)/2 — axi/r. For H = (|5^ + r^)/2 one can take F = {H; IIi,Il2, . ■ ■ , Hn-i, 
Pi2, Pi3, ■ ■ ■ , Pin), where Hi = \(pf + xf). For H = f{P^) one can take 
F = (P^;p^, r, L). In all cases except the first one, we have k = 1. 

Since x ■ p — {p^r^ — r'^p^)/A — n/2, this proposition implies in particular 
that [x ■ p,P] = 0. Note also the relation — r'^p^ — [x ■ p)'^ — (n — 2)x ■ p, 
which can for instance be easily verified using proposition 3.2 of [T]. 

Proof. Let the operator H of class Ok have the form H = f {p^ , P"^ , gi{r) , 
. . . , gi{r)), where the functions f,gi, . . . ,gi have the properties specified in the 
proposition. Then the relation [i? , P] — follows from (12.131) and (j2.14p . 

Viceversa, let H be an operator of class O such that [H , P] = 0. Then 
lemma 3.24 of [T] implies that {MH, P} — 0, where H and P are the symbols of 
H and P respectively. Using lemma 1^?^ we thus obtain that MH — g{p^, r,P^), 
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where g is an arbitrary function of three variables. Furthermore, since and 
are both homogeneous polynomials of order 2 in p, taking into account the 
definition of main part we obtain that g is a homogeneous polynomial in the 
two variables {p^, P^), whose coefficients are arbitrary functions of r defined for 
all r > 0. 

The proof of the remaining statements is similar to the proof of the cor- 
responding statements of proposition 12.31 Let us consider, in particular, the 
hamiltonian H = p'^ /2 — a/r oi the quantum Kepler system [TS]. Using (|2.13l) . 
together with proposition 2.1 and lemma 2.2 of [5], it is easy to verify that 
[H,A^] = 0, where 

n 

= ^ Pij o pj - a— , i = l,...,n. 
i=i ^ 

We have used above the symbol o to denote symmetrized products, as in [5]. 
Only one component of A is quasi-independent of the set {H,P^,L). Note that 
we have in this case 

= 2H 

Finally, it is easy to check that the main parts of the operators of the sets 
considered in the last part of the proposition are functionally independent. □ 

The commutation relations between operators, from which the integrability 
of the considered sets of operators has been established, have been derived 
exploiting the quadratic dependence of classical momenta Pij on the canonical 
variables {x,p). Let us now present an alternative proof of these relations, which 
is only based on the linear dependence of these momenta on impulses p. We 
shall consider the quantization of an arbitrary vector field on configuration space 
K, more exactly, the quantization of the hamiltonian function on T*K which 
corresponds to the lifting of this field on T*K. These considerations are useful 
for the investigation of any linear operator which is invariant with respect to 
the phase flows of such vector fields on K, independently of the assumption that 
these vector fields be linear. 

Let P = (Pi, . . . , P;) be a set of functions on the symplectic manifold M — 
T*K, which are linear with respect to p: 

P, = v^{x) + {p,v,{x)), (2.15) 

Pi : T*K — > M. Such functions, in local coordinates (x,p) induced by local 
coordinates x on K, have the form Pi = Pi{x,p) — v^{x)+J2k=i ''^i{^)Pk- Let us 
suppose that the linear combinations of these functions with constant coefficients 
form a Lie algebra g with Poisson brackets in the role of commutators, and that 
the functions of the set P form a basis of this algebra. Let us consider the set 
of operators P = (Pi, . . . ,Pi) obtained by standard quantization from the set 
of functions P, i.e., 

A :-t'."(x)+^i;f(x)— , i = l,...,l. (2.16) 
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Then the Hnear combinations of these operators form also a Lie algebra with 
respect to the usual commutator of linear operators. Moreover, let us consider 
the linear map, from the original Lie algebra of functions on M to the Lie algebra 
of operators, which is defined by the correspondence of sets P P, obtained 
by standard quantization. This map is a isomorphism between these two Lie 
algebras, i.e., the linear map preserves commutators. 

This fact is an obvious consequence of the following more general proposition. 
Let us consider the Lie algebra V = Vect k{K x E) of all vector fields defined on 
the direct product KxM.3 {x,u), which do not depend on u. The commutator in 
this algebra is the usual Lie bracket of vector fields. In local coordinates (x, u), 
X = (xi, . . . ,x„), a vector field V Q V on the (n + l)-dimensional manifold 
if X M has the form x — v{x), u — v^{x). Let us consider also the algebra F 
of all functions P on the symplectic manifold M = T*K which are linear with 
respect to the impulse p, i.e., functions of the form (j2.15p . Let us also consider 
the Lie algebra O of all linear nonhomogeneous differential operators on i.e., 
operators of the form (j2.16p . 

Proposition 2.6. There are canonical isomorphisms between these three Lie al- 
gebras V, J^, C In local coordinates on K , the coefficients Vq{x)^ vi(x)^ . . . , Vn{x) 
defining the elements of these algebras are conserved under these isomorphisms. 

Let V and W be two vector fields of class V = Vccti<-(i4r x R). Let us indicate 
the functions and operators, associated with these fields, as Py, Pw and Py, Pw 
correspondingly. Then the conservation of commutators of these algebras under 
the considered isomorphisms can be written in the form 

{Pv, Pw} = P[v,w] , [Pv, Pw] = P[V,W] ■ 

The hamiltonian vector field Xp, defined by the hamiltonian function P G 
T , can be lowered by natural projection n : T*K — > K. This means that 
Tr^:{Xp{m)) G Txl^ does not depend on the choice of the point m S it~^(x), 
where Xp{m) is the vector field Xp at point m, tt, : TM — >■ TI'C is the deriva- 
tive of the map n, and T^K is the tangent space to K at point x. Suppose that 
P = Py, where V € V. Since the elements of V do not depend on u, the vec- 
tor field V can be lowered onto configuration space K via the natural projection 
K X R„ — )• K . These two vector fields, obtained by projection on K from Py 
and V respectively, are coincident. 

Proof. Both statements of this proposition, about the correspondence of com- 
mutators of the three Lie algebras and about the coincidence of the projections 
on K of the two vector fields, can be easily checked by direct computation 
in local coordinates. In general, these statements are reformulations of simple 
well-known facts. □ 

Since classical momenta Pij{p,x) = XiPj — XjPi are linearly dependent on 
classical impulses p, we can use the above proposition to deduce the commu- 
tation relations (|2.12p for the operators Pij from the corresponding classical 
relations (|2.5p . 

Relations (|2.5p show that the set of momenta Pij, 1 < i < j < n, is a 
basis of the Lie algebra g = so(n), which is isomorphic to the Lie algebra of 
all skew-symmetric matrices. A natural basis in this Lie algebra is formed by 
matrices D^^ , ^ < i < j < n, defined by formula (|2.2p . The correspondence 
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£)u K- > is extended to linear combinations of matrices D'^^ and functions 
Pij as an isomorphism of Lie algebras. The basis P = (-Py , 1 < i < j < n) 
of Lie algebra g induces a dual set of coordinates on the co-algebra g*. It 
is well-known that the function : g* — > M, P^ := J2i<i<j<n ^ij^ 
invariant of the co-adjoint representation of SO{n) on the co-algebra q* , where 
g = so{n). From corollary 4.4 of [5], it follows that [(P^)p"\P] = 0, where 

(P^)p™ denotes the symmetrization with respect to P of polynomial P^. But 
obviously (P2)J'" = Ei<,<j<„^j = P^- We thus conclude that [P^P] = 0. 

Let us prove now that [t'^jP] = [p^,P] = 0. Let us consider the set 
B = {l,x,p,P) of Z := n(n-l)/2-f2n + l functions on M^" = T*K. From Pois- 
son brackets relations (j2.5p . (j2.7|) and (|2.8I) it follows that this set of functions 
is a basis in a /-dimensional Lie algebra. The functions of set B are linear non- 
homogeneous functions of p. Therefore proposition 12.61 implies that analogous 
commutation relations hold for the operators B — (1, x,p, P) obtained from the 
functions of set B by standard quantization. We have in particular 

[xi,Pjk\ = 5i]Xk - kkXj , (2.17) 
[Pi,Pjk] = Si-jPk - SikPj . (2.18) 

Then relations (j2.13p can be derived from (j2.9p using proposition 4.2, case b, of 

The following proposition is the quantum equivalent of proposition 12.41 

Proposition 2.7. For any n > 2 and for any z — 1, ... ,n — 1, it is possible 
to construct in a recursive manner sets Zn^z o,nd Ln.z of polynomial functions 
of degree < 2 in the variables P„ :— (Py , 1 < i < j < n), with the following 
properties: 

1. Zn_z contains z elements, 

2. Ln_z contains 2{n — z — 1) elements, all of degree 1 in P„, 

3. the set Yln.z '■— {Zn.z, Ln,z) is quasi-independent, 

4- [^n,z, ^n.z] = 0. 

Proof. Let Zn_z and L„^z be the polynomials obtained from the classical ones 
Zn.z and Ln,z of proposition 12.41 by simply replacing their arguments Pn with 
the quantized momenta P„. In this case symmetrization is unnecessary, since in 
these polynomials all monomials of degree 2 are squares of elements of P„ . Since 
all elements of the set Iln,z = {Zn.z, Ln,z) are homogeneous polynomials in p, 
they coincide with the symbol of the main parts of the corresponding elements 
of tln,z- Therefore the quasi- independence of II„^z follows from the functional 
independence of the classical set Iin,z- 

Point 4 can be proved just by repeating the proof of proposition 12.41 When 
the set of indexes N,-a :— (l,...,m) is split into two disjoint subsets Ii and 
I2, consider in fact the two sets of operators P'-'^^ C Pm, k — 1,2, which are 
the standard quantization of the sets of momenta P*^*^-* . From the isomorphism 
between the two Lie algebras generated by the sets P„ and P„ respectively, it 
follows that [P(1),P(2)] ^ 0. Moreover, since P^ is a Casimir function for the 
co-algebra so(m)*, from coroUary 4.4 of |2 it follows that [P^^, Z„j^^^j^ (P^^^)] = 
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[P^, in^^zj. (P*^*^^)] = for fc = 1,2. Hence point 4 is obtained by induction on 
rt, as in the classical case. □ 

Let H be an arbitrary operator of class O owK — ]R"\{0}, such that [H, P] — 
0. From proposition 12.51 it follows that MH — g{p^ ,r, P^), where the function 
g is a homogeneous polynomial in the two variables (p^jP^), with coefficients 
dependent on r. If the polynomial g satisfies the condition dg/d{p'^,r) ^ 
almost everywhere, from proposition 12.71 it follows that the sets of operators 
Fn,z ■= {H , Zn.z'i Ln,z) are quasi-integrable sets, with subset of central elements 
{H, Zn,z), for all 2; = 1, . . . , n — 1. We have ^Fn,z = 2n — k and k = z + 1. Hence 
the number k of central elements can take all values from k = 2 to k = n. 



3 Free rotation of an n-dimensional rigid body 

The configuration space K of this system is the group SO{n) of orthogonal 
transformations of the 7i-dimensional euclidean space K = SO{n). Its phase 
space is Af = T*K = T*SO{n). We shall present a detailed analysis of the 
integrability of the classical system, which distinguishes itself from other already 
existing investigations [TBI IIZl HH] , in the fact that we apply here the concept 
of noncommutative integrability (see definition 3.16 of [1 ), and analyze the 
dependence of the number k of central integrals on the properties of the set 
of generalized moments of inertia, more precisely on the presence in this set 
of subsets consisting of moments equal to each other. The results will then be 
applied to the investigation of the integrability of the corresponding quantum 
system. We begin however with a preliminary subsection on some properties 
of the group SO{n), most of which are more or less well-known, but which 
shall here be derived in a form and with a notation convenient for our present 
purposes. 

3.1 Properties of left- and right-invariant vector fields on 

SO{n) 

Any arbitrary Lie group G acts on itself by left and right shifts: for each element 
g (z G these are defined by the diffeomorphisms 

Lg : G G, Lgh — gh , 
Rg-.G-^G, Rgh:^hg. 

Let us denote by TgG the tangent space to G at point g. Then the Lie algebra 
g associated with G can be identified with the tangent space g = TgG at the 
neutral element e of the group G. Let {Lg)^ and (Pg)* respectively denote the 
derivatives of the maps Lg and Rg at e. Each element a S g defines two vector 
fields and on G. At any point g & G these vector fields are respectively 
defined as V^^{g) := {Lg)^a € TgG and Vj^{g) (Pg)*a e TgG. We have 

(L,),K'^(.9) = {UULgU = iUgU = V^(hg) , 

{Rh)*V^{g) - {RhURg)*a = {RgnU = V^{gh) . 

These relations mean that the vector field (respectively V^) is invariant 
with respect to the action of Lie group G on itself by left (respectively right) 
shifts. 
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Definition 3.1. For the above reason, the fields and are respectively 
called left-invariant and right-invariant vector field on G associated with the 
element a £ g. 

The Lie brackets of left-invariant vector fields respect the structure of the 
Lie algebra g: [V/", V/"] — Vj^^j, where a, 6 e g and [a, b] is their commutator in 
Q. An analogous result, although with a reversed sign, is true for right-invariant 
vector fields: [V^,V^] = -V^a.by We have also [V^ ,V^] = 0. 

The Lie group G = SO{n) can be identified with the group of orthogonal 
matrices of size n x n, and its associated Lie algebra g — so{n) with the Lie 
algebra of skew-symmetric matrices of the same size. Then the vector field 
(respectively V^) corresponding to the skew-symmetric matrix A is defined by 
the system of differential equations X = XA (respectively X = AX)^ where 
X is an orthogonal matrix. The commutator of the Lie algebra g — so{n) is 
given by the usual commutator of matrices. We can therefore rewrite the above 
formulas for the Lie brackets between left- and right-invariant vector fields as 

[Vl,Vi]^~V^ls], [V^i',V#]=0. (3.1) 

We can also consider the classical impulses P^,P^ : T*M — > K, with AI — 
SO{n), which are associated according to proposition 12.61 with the vector fields 
Vj' and Vf, for any A e g: 

Pi(m) ^ {p,ViiX)) , P^im) ^ (X)) , (3.2) 

where p = m G T^G. These impulses form a Lie algebra with respect to 
Poisson brackets, which is isomorphic to the Lie algebra of the corresponding 
vector fields. Therefore, from relations (j3.1[) we derive 

{Pk,Pk} = P[lB], {Pa,pE}^-P[Ib], {Pi,P|} = 0. (3.3) 

Let us consider the vector fields V^^, V;f, which are associated with the 
matrices D*^ defined by formula (|2.2p . We recall that the set (-D*-' , ^ < i < j < 
n) forms a basis of the Lie algebra g = so{n) = T^G. This implies that V^{g) 
and Vjj{g) are two bases of linear space TgG, for any g G G. From (12. 3p and 
p.ip we obtain that the Lie brackets between these vector fields have the form 

-<5.,.V^t - S,kV^ + 6,kV^i + 5,hV,1 , (3.4) 
S^hV^1 + 5,uVi - kkV^r. - SjhVi , (3.5) 
0. (3.6) 

These vector fields are associated with the classical impulses 

P^(m) = (p, V^{X)) , /^f (m) = (p, V^{X)) . (3.7) 

According to (|3.3I) , the Poisson brackets between these functions have the same 
form as the Lie brackets (j3.4D - (j3.6p between the corresponding vector fields: 

-S^HPji - 5,kPk + 5^kP^H + 5jHP^tk , (3.8) 
5^hP^k + S.kPi - S^kP^h - , (3.9) 

0. (3.10) 



{Pij^Phk} = 
{^if I Phk\ — 
{Pij^ Phk} = 
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In the following we shall often indicate with = P^{m) and — P^ijn) 
the two skew-symmetric matrices having components P^{m) and P^{m) respec- 
tively, with 1 < i,j < n. Since for any A G we have A = X]i<j ^ijD^'' ^ from 
(1321) we obtain P^ = E.<, A^.P^j = ~\ E.,, Pt,A,^ = {P^A). Using the 

first of relations p.3p we then obtain 

{Pk,PL} = -iTr(P^[Ai?'']) - -iTr(i?'''^[P^A]) = [P\AU, (3.11) 

where the last member represents the element at row h and column k of the 
commutator of the two matrices P^ and A. In a similar way we obtain 

{P^.P^,}^-[P^,AU. (3.12) 

Remark 3.1. Since the N vectors l^^(X), 1 < « < j < form a linear basis of 
TxG for any X £ G = SO{n), we see from formula (|3.7p that the correspondence 
p i-> P^{m) is an isomorphism between the linear spaces T^G and so{n). It 
follows that, for a generic m G M, the matrix P^{m) is a "typical" n x n 
skew-symmetric matrix. The same fact is obviously true also for P^{m). 

Let us represent the generic element g of the group G = SO(n) as an or- 
thogonal n X n real matrix X . If we indicate as X the transposed of the matrix 
X, so that Xp^ := X^js, the orthogonality of X implies X~^ = X. For any 
skew-symmetric matrix A — —A, at any point X G SO{n) the tangent vec- 
tor X = AX can also be written as X — XB, where B = X~^AX is also a 
skew-symmetric matrix. Using this fact, it is easy to see that the basic left- 
and right-invariant vector fields introduced above are connected to each other 
by the relations 

n 

V^{X)= X,hX,^VtAX). 

h,k=l 

or in matrix notation 

V'^iX) = XV^{X)X . (3.13) 
Similarly, we have for the corresponding impulses 

n 

P"i^)= E X,hX,kPL{m), (3.14) 

h,k=l 

or equivalently 

P'^im) ^ XP^{m)X , (3.15) 

where m G T^G. 

It is sometimes useful to indicate with and P^ also the two sets of 
N = n{n — l)/2 functions on T*G defined by formulas (|3.7p : 

P^ = (f^,l<«<J<«), P'' = {P,'fA<i<J<n). (3.16) 

In this way the elements of the sets P^ and P^, evaluated at a point m G M, 
just coincide with the independent elements of the two skew-symmetric matrices 
P^{m) and P^{m) respectively, which were introduced before formula (j3.11l) . In 
the following, it should be clear from the context whether P^ (or P^) indicates 
a skew-symmetric matrix or the corresponding set of N independent elements. 
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Let F = {Fi, . . . , Fr) be a set of r real functions Fi : M — > M, i = 1, . . . , r, 
on a manifold M. As usual, we denote as rankF at m e M the dimension 
of the image F^,{TmM) of tangent space T,nM with respect to the derivative 
F^ of the map F : M ^ W defined by the set F. Let dF = {dFi, ...,dFr) 
denote the set of the differentials of the elements of F at the point m. It is well 
known that ranki^ = dim Span di^, where SpandF denotes the linear subspace 
of T^M spanned by the elements of dF. 

Let M = M"^^ be a symplectic manifold. Let us denote with H the antisym- 
metric bilinear functional H : T*^M x T^M — M such that the Poisson bracket of 
any two functions /, /i : M — K is given at m by the relation {/, h} = n((i/, dh). 
For any linear subspace L C T^M, let us denote with the subspace skew- 
orthogonal to L with respect to H, i.e., {m G T^M : Il{u, w) = Vw S L}. 
Since H is nondegenerate, we have dimi -I- dimL^ = 2N. From these facts we 
deduce the following 

Lemma 3.1. Let F be a set of functions on a symplectic manifold M"^^ . Then 

dim (SpandF)^ = 2N - ranki^ . (3.17) 

Proposition 3.2. Let us consider the two set of functions and defined 
by formulas iS.lb]] and \3. 7p . On all M ~ T*G, where G — SO{n), we have 

rankP^ = rankP^ . (3.18) 

In addition, the two subspaces SpandP^ and SpandP^ of T^M are related to 
each other by the equalities 

(SpandP^)^ = SpandP^ , (SpandP^)^ = SpandP^ . (3.19) 

Proof. Let us fix a system of coordinates x in a neighborhood of g € G, and 
let {x,p) be the local coordinates induced from x on T*G. According to p.7p . 
we have V^^ — dP^ / dp. Since at any point 5 of G the N vector fields , 1 < 
i < j < n, are linearly independent, we have that the differentials dpP^ £ TgG 
with respect to variables p of the N functions P^j are linearly independent. This 
implies in particular the independence of their differentials dP^^ G T^M with 
respect to variables (a;,p), so that rankP^ = TV. We obtain in a similar way 
that rankP^ = TV. Equalities p.lSp arc thus proved. 

Formula p.lOp implies SpandP^ C (SpandP^)^. Furthermore, formula 
p.lSp is equivalent to dim Span dP^ = dim Span dP-'^ = N . Therefore, applying 
lemma [33] to the set P^ we obtain dim(Span(iP'^)^ = N = dim Span dP^. 
These facts imply the former of equalities p.l9p . The latter is obtained in a 
similar way. □ 

Let B denote the set B := {P^,P^) of 2N = n{n - 1) functions on T*G. 

Corollary 3.3. On all M = T*G we have 

(Span dP)^ = Span dP^ n Span dP^ . (3.20) 

Proof. Recalling equalities p.l9p . from the relation B = P^ U P^ we obtain 
(SpandP)^ = (Span dP^)^ n (Span dP^)^ = Span dP^ n Span dP^. □ 
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In order to establish the hnear dimension of subspace p.20p . it is first nec- 
essary to recall a basic property of skew-symmetric matrices. 

Lemma 3.4. Let us consider a skew- symmetric operator in an n- dimensional 
euclidean space L. Then there is a cartesian system of coordinates in which the 
operator is defined by a matrix which has the "normal block- diagonal" form: 



A = 



/ ai 
-ai 












^ 








Q!2 
-0(2 
































as 
-a, 
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) 



(3.21) 



where s = [^] is the integer part of ^, at > for k ~ 1, . . . , s, and the last 
vanishing row and column are present only for odd n. Let us suppose that all the 
eigenvalues ±iai, . . . , ±ias, (0) of A (where i — ) are pairwise different. 
(This means that au > 0, at ^ ak^hjk — l,...,s, h ^ k). Then a skew- 
symmetric matrix B commutes with A, i.e., [A, B] — 0, if and only if B has the 
same block- diagonal form, without any condition on its eigenvalues. 

An equivalent intrinsic formulation is the following. For any skew- symmetric 
operator A with pairwise different eigenvalues on the euclidean space L, there ex- 
ists a unique decomposition i = L^SL^® ' ' '©^^©-^0 '-'f ^ ^^^^ ^ 2- dimensional 
invariant subspaces Lf and (for odd n) a 1-dimensional subspace Lq. Any skew- 
symmetric operator B commutes with A if and only if all subspaces if 
of this decomposition are invariant under the action of B . 



id Ln 



We shall denote with gtyp C g = so{n) the set of all skew-symmetric matrices 
whose eigenvalues are pairwise different. Obviously, "almost all" elements of g 
also belongs to fltyp- More exactly, g\gtyp is a closed subset of g having vanishing 
measure. It follows that, if some statement is true in gtyp, then it is true (at 
least) almost everywhere in g. 

Let us fix an element a of the Lie algebra g, and consider the linear operator 
ada : g — 7> g in this algebra, ada : 6 i->- [a, &]. The kernel of ada is the subalgebra 
of all elements of g which commute with a, i.e., Kerada = {b G g : [a,b] — 0}. 



Corollary 3.5. If a £ gtyp, the 



dim Ker ada 



(3.22) 



Let ada(g) C g denote the image of g with respect to the operator ada- Then 



dim ada (g) ~ N 



(3.23) 



Let a be represented by a matrix A having the normal block-diagonal form h3.21]) . 
with ak > 0, ah ^ ctk"^ h, k ~ 1, . . . , s, h ^ k. Then ada(g) is the linear space 
of the matrices B d g such that B12 = -B34 = • • • = B2s-i.2s = 0, with s = [n/2\, 
i.e., 

ad^(g) = {Beg: B2^-l,2^ = Vi = 1, . . . , [n/2]} . (3.24) 
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Proof. Equality (I3.22p follows from lemma [5^ and from the observation that 
the linear space of the matrices of the form p.2ip has dimension s — [2^] . 
Equality (j?^ follows directly from (jX^ . since dimg = N = n{n - l)/2. 
A direct computation shows that, for any C € q, one has -621-1.21 = Vi = 
1, . . . , [^] , where B :— [A, C]. Therefore the linear space on the right-hand side 
of (|3.24l) contains ad^(g). Furthermore, it is obvious that the right-hand side 
of (I3.24[) has linear dimension N — [^] . Taking into account (|3.23p . these facts 
imply ((XMl) . □ 

For a given point m £ M = T*SO{n), let us consider the two skew- 
symmetric matrices P^{m) and P^{m) whose elements are defined by formula 
(13. 7|) . Let us denote with Mtyp C M the set of all points m £ M such that 
P^{m) g gtyp. We see from formula (I3.15P that P^{m) S gtyp if and only if 
P^{m) S 0typ- We also know (see remark [3TT|) that, for a generic m € M, the 
matrix P^{m) is a typical nxn skew-symmetric matrix. Hence, the eigenvalues 
of P^{m) are pairwise different for almost all m G M . It means that almost all 
elements of M belong to Mtyp, or equivalently that M \ Mtyp is a closed subset 
of M having vanishing measure. 

Let (7 — (7 (to) be the linear dimension of the kernel of the linear operator 
adpi : 

a := dimKeradpt . (3.25) 

In other words, a is the linear dimension of the subalgebra of all matrices A £ 
so{n) such that [P^,A] = 0. 

Lemma 3.6. For all to e AI , where M — T*SO{n), we have 

dimKeradpK = dimKeradpt = a . (3.26) 
In addition, for all to G Mtyp (hence, for almost all to G M ) we have 

Proof. According to formula (j3.15L VA e g = so{n) we have 

[P^,A] = [X-^P^X,A] ^ X-^[P^\kdx{A)]X , 

where Adx(^) = XAX~^ is the image of A with respect to adjoint action 
Adx : — >■ QoiX ong. Hence, A G Keradpt if and only if Adx (^) S Keradpi?.. 
This implies p. 261) . since the adjoint action Adx is an algebra automorphism 
of g. Equality (13:27)) follows from formula ((3?22|) of corollary [331 □ 

Proposition 3.7. Let us consider the set B = (P^, P^) of2N functions defined 
by formula i3.7\ l. At any point m £ M = T*G, where G — SO(n), we have 

(SpandB)^ = |a e T*,M : a = ^ AjdP.^j , A^-A, [P^,A] = o| 

i<j 

= {a e T^M : a = ^ A.jdP,^ ,A^~A, [P^ , A] = o| (3.28) 

and 

ranks = 2A^- cr, (3.29) 
where a is defined by formula iS. 25]} . 
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Proof. Equalities ((XTOl) and l!^^ imply that 



(SpandB)^ = SpandP^ n (SpandP^ 



Therefore (Spandi?)^ is the set of all those elements a E SpandP^ such that 
n(a, dP^j.) = V/i, fc with 1 < h < k < n. According to (IXTS)) . dP^ is a set 
of linearly independent elements of T^G. Hence any covector a € SpandP^ 
can be expressed as a = J2i<j ^ijdPij, where Aij are elements of a univocally 
determined skew-symmetric matrix A. We have 

U{a, dPt,) - A,MdPtj,dPtk) = {Pa, PL} = [P\AU , (3.30) 

where for the last equality use has been made of formula (|3.1ip . Therefore 
a G {Span dP^)^ if and only if [P-'^,j4] = 0. This implies the first equality of 
formula (j3.28l) : the second one can be obtained in a similar way. 

Formula (j3.28l) shows that (Spandi?)"^ is in one-to-one correspondence with 
the subalgebra of matrices A € so{n) such that [P^, A] = 0. Hence 

dim(Span(iP)^ =dimKeradpj^ = f , (3.31) 

so that (I3.29P is obtained by applying lemma IXTl to the set P. □ 

From formulas (I3.29P and p.27p we immediately obtain the following: 

Corollary 3.8. At all points of Mtyp (hence, almost everywhere on M ~ T*G) 
we have 



rankP = 27V 



(3.32) 



The previous corollary, together with proposition l2.11 implies that there exist 
at most [^] functionally independent functions on M which are in involutions 
with the whole set B. In order to explicitly construct a set of [^] such functions, 
it is useful to extend to Lie algebras the notion of Casimir function which has 
already been introduced for Lie co-algebras in [2]. 

Definition 3.2. A Casimir Junction on a Lie algebra q is an invariant of the 
adjoint action of the local Lie group G, i.e., a function which is constant on the 
orbits of this action. In other words, a Casimir function i^T : g R is a first 
integral which is common to all differential equations in g of type b = [a, 6], 
where & G g and a is a fixed element of algebra g. 

Let us consider the case g = so{n), i.e., g is the algebra of skew-symmetric 
matrices. Fix A E Q and consider the characteristic polynomial det(A — AP), 
where E is the identity matrix. The relation A = —A implies det(A -|- \E) = 
det(A -I- \E) = (— l)"(y4 — XE), so that the characteristic polynomial has the 
form det(AP - A) = A" + Ci{A)X"~^ + C2{A)X"-* + The coefficients 
Ci{A), . . . ,Cs{A) of this polynomial, where s = [^] , are clearly polynomial 
functions of the N = ri(n — l)/2 independent elements Aij, 1 < i < j < n, of the 
skew-symmetric matrix A. Such elements are the coefficients of A with respect 
to the basis D'^ of g given by formula (gSl). We will caU C = (Ci, . . . , C,) the 
standard set of Casimir functions on so(n). 
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Lemma 3.9. The functions of the standard set C of Casimir functions on 
Q = so(n) are really Casimir functions in the sense of definition \3.'A Moreover, 
the functions of set C form a basis in the space of all Casimir functions in 
the following "functional" sense. The differentials of the functions of set C are 
linearly independent at all points o/gtyp' therefore, they are linearly independent 
at almost all points ofg. Moreover, every Casimir function K on q can be locally 
expressed as a function of the elements of this set: K = K{C). 

Proof. These are actually well-known facts, but we give the proof for complete- 
ness. The adjoint action of the group G = SO{n) on the Lie algebra g = so(n) 
takes in the matrix representation the form A i— Adx(^) = XAX~^, where 
X e G, A e g. Therefore det(yl- A£;) = det{Adx{A)~ XE), i.e., the polynomial 
det(^ — XE) in the variable A is invariant under the adjoint action of G. This 
means that the coefficients of this polynomial, i.e., the functions of set C, are 
also invariant, and are thus Casimir functions in the sense of definition 13.21 

Let us now examine the functional independence of the s functions of set C, 
where s = [§] , at a given point B e q. To this end, we shall evaluate the rank 
of the 3 X N matrix 

, 1 <h < s, I < i < j < n) , 
B=B J 

whose rows are labelled by the index h and the columns by the double index if. 
In the above formula, we have denoted with dCh{B) / dBij the partial derivatives 
of the Casimir function Cu with respect to the independent elements Bij of the 
skew symmetric matrix B. We know from lemma 13.41 that there exists X — 
X{B) £ G such that AAx{B) = A, where A is a skew-symmetric matrix in the 
normal form p.2ip . such that > for fc = 1, . . . , s, and ±mi, . . . , iia^, (0) 
are the eigenvalues of B. Since C is a set of Casimir functions, we have that 
G{B) = G{A). Taking into account that Adx is an invertiblc linear operator in 
0, it follows that rank J(B) = rank J(A). 

Let P{B, X) := det(A£' — B) denote the characteristic polynomial for the 
matrix S G g. For any A G M it is easy to see that 

= Vi, j = 1, n such that j > i + 1 . (3.33) 

B=A 

Recalling the definition of the set C, equality (j3.33p implies that all the columns 
of matrix J {A), corresponding to indexes ij with j > i + 1, are zero. We can 
thus write rank J(A) — rank J(A), where we have introduced the s x s matrix 

, 1< h< s, l<i<s) . (3.34) 

B=A / 

It is easy to see that, for a matrix A of the form (j3.21L the character- 
istic polynomial has the form det(A£' - A) = (A)ni=i('^^ + = A" + 
^ - af A"~^ -I- J2i<j '^j'^j -^""^ + Therefore the standard set of Casimir func- 
tions can be written as G{A) = D{/3{a)), where a = {ai, . . . , as), [3(a) — 
(/3i(a), . . . , A(a)) (a?, . . . , a^), D{(3) = (7^i(/3), . . . , DM), and 

DhiP):= AiA2---ft. for h^l,...,s. 

il<il<---<ih 



J = J{B) 



dCh{B) 



dP{B,X) 
dB,, 



J.J,,,,. 
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Note that in particular Di{/3) — J2i=iPi ^^'^ Ds{(3) — Pil32---(3s- Since 
A2i-i,2i = a^, we have 



dCh{B) 



2j-l,2i 



= Za, 



B=A 



Therefore, recahing p.34L we have det J(A) — 2^aia2 ■ ■ ■ asKs{f3{a)), where 
Ks{f3) is the determinant of the s x s jacobian matrix It is easy to 

see that Ks{P) is a symmetric polynomial of order s(s — l)/2 in the variables 
/3, which vanishes whenever = /3j for some i ^ j. Therefore we can write 
KsijS) = CsWj^^jiPi — f3j). In order to evaluate the constant coefficient Cs, 
we note that, when fig = 0, one has dDg/dPi = for i = 1, . . . , s — 1, and 
dD,/dps = Pih ■ ■ ■ Ps-i- It easily follows that A',(/3) = /3i/?2 • • • (is-iKs-i{fi') 
for = 0, where /3' = . . . ,/3s_i). Using this fact, it can be easily shown, 
by induction with respect to s, that = IVs e N. We thus conclude that 

det J{A) = 2"aia2 • • • J|(a,^ - a^j) . (3.35) 

When all eigenvalues of matrix B are pairwise different, we have > 0, 
ah ^ ak.y h,k — 1, . . . ,s, h ^ k. In such cases we see therefore from p.35p 
that det J(A) ^ 0, so that rank J(i3) = rank J(^) — s. This means that the 
differentials of the s Casimir functions of set C are linearly independent at B. 

In order to complete the proof of the lemma, we note that, according to 
corollarv l3.5[ at almost each point B G g the set of vectors [A, B], where A varies 
on all 0, forms a subspace of codimension s = [^] in g. Since the differential of 
any Casimir function must be zero when acting on this subspace, we deduce that 
the differentials of any r Casimir functions, where r > s, are linearly dependent 
at any point B d g. Therefore the differential of any Casimir function K at any 
point of algebra so{n) is a linear combination of the differentials of the functions 
of the standard set, i.e., we have locally K — K{C). □ 

Note that the set of all Casimir functions on g defines the orbits O of the 
adjoint representation of the corresponding local Lie group G by their common 
level surfaces: C~^(c) — O, where C~^(c) is the pre-image of a point c € W, 
and s = [^] . More exactly, this is true in the domain of regularity of the map 
C : g W, defined by the set C. From this it is easy to obtain that the typical 
orbit O has codimension s in g, i.e., dimO = dimg — s. 

For g — so{n), we can associate with every element & G g* a unique skew 
symmetric matrix B such that {A,b) — ^Ti (AB) = J2i<j AijBij VA G g. 
Clearly, the elements Bij, I < i < j < n, oi matrix B are just the coefficients 
of b with respect to the dual basis of the basis D^^ of so{n) defined by formula 
(|2.2p . According to this one-to-one correspondence, in the following we will 
often identify g* with the co-algebra of skew-symmetric matrices, and we shall 
denote with g*yp the set of all elements of g* whose eigenvalues are pairwise 
different. Obviously, almost all elements of g* also belongs to gjyp. It is easy 
to see that, in the matrix representation, the co-adjoint action of the group 
G — SO{n) on an element A G g* has the same form as the adjoint action on 
the corresponding element of g, i.e., A i— >■ XAX~^, with X G G. It follows that 
the standard set of Casimir functions on the Lie algebra so(n), introduced in 
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definition 13.21 also represents a standard set of Casimir functions on the space 
of skew-symmetric matrices, when the latter is identified with the co-algebra 
so{n)* . This can be formally stated as follows. 

Definition 3.3. Let g* be the dual co-algebra of g = so{n). Fix a skew- 
symmetric matrix A d q* and consider the characteristic polynomial det{XE — 
A) = A"-HCi(A)A"-2-hC2(A)A"-4-F-- - . The coefBcients Ci(A),... , C^(A) of 
this polynomial, where s = [^] is the integer part of ^, are polynomial functions 
of the set a of the matrix elements of A. We will call C — (Ci, . . . , Cs) the 
standard set of Casimir functions on so{n)* . 

Lemma 3.10. The functions of the standard set C of Casimir functions on 
the coalgebra q* , where g = so{n), are really Casimir functions in the sense of 
the definition given in J^. Moreover, the functions of set C form a basis in the 
space of all Casimir functions in the following "functional" sense. The differ- 
entials of the functions of set C are linearly independent at all points of g^y^ . 
Therefore, they are linearly independent at almost each point of q* . Moreover, 
every Casimir function K on q* can be locally expressed as a function of the 
elements of this set: K — K{G). 

Proposition 3.11. Let C be the standard set of Casimir functions on the Lie 
algebra so{n). Then C{P^) = C{P^) on all symplectic manifold M = T*G. 

Proof. From relation p.isp . taking into account that X ~ X^^VX E SO{n), 
we obtain that det(A£' - F^) = dct{\E - P^) on aU M for any A e R. The 
thesis then follows from definition 13.91 of the set C. □ 



In the following we shah denote with C the set C := C{P^) = C(P^) of 
s = [^] real functions on M . 

Proposition 3.12. On all M ~ T*G, where G — SO{n), we have 

{C,B}m-0. (3.36) 

Moreover, on all Aftyp (hence, almost everywhere on M ) we have 

rank (7 



n' 
2. 

z 



(3.37) 

Span dC = (Span dB)^ . (3.38) 

Proof Since C = C{P^), formula ((XTU)) implies {C,P^} = 0. On the other 
hand, according to proposition 13.111 we also have G = G{P^), so that using 
again (j3.10l) we obtain {C, P^} = 0. Equality (j3.36l) is thus proved. 

Since both g and g* have been identified with the set of n x n skew-symmetric 
matrices, we have Mtyp = {m G M : P^{m) e fltyp}- Lemma [3.91 implies that 
rankC = [^] at all points of gtyp. Moreover, according to (I3.18p . P^ : M g 
is a regular map. Therefore, it follows from the definition C ~ G o P^ that 
equality p.37p holds at all points m € -Mtyp- 

Formula (j3.36p implies that SpandC C (SpandS)^. Moreover, from p.37L 
((33T|) . and ((3?27l) . we obtain that dim Span = [f] = dim(SpandB)^ at aU 
points of Aftyp. These relations imply (j3.38p . □ 
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Taking into account p.l9p and p.20p . equality p.38p also implies 



SpandC = Span (iP^ n SpanrfP^ = SpandP-^ n (SpandP^)^ . (3.39) 

Remark 3.2. For any function f : M ^ M, the condition {/, P} = is equivalent 
to df £ (SpandP)^. According to formula (j3.38p . the latter condition implies 
that locally / = g{C), where g : R["/^1 — )■ M is an arbitrary function. Recalling 
the definition of the set C, this means that / = g[C{P^)) = c{P^), where C is 
the standard set of Casimir functions, and therefore c := g o C is also a Casimir 
function. Taking into account proposition 13. Ill we also have / = c{P^). We 
conclude that the functions on M which arc in involution with the whole set B 
are all and only the Casimir functions of the set P^ (or equivalently P^). 

3.2 Free rotation of a classical rigid body 

Let us consider the system with hamiltonian function of the form 



where J is an operator in the A^-dimensional linear space of impulses P^, 
with N ^ n{n- l)/2. The map P^ : T*G ^ = fl is linear with respect 
to impulses p, see formula p.7p . Therefore the function H = H{P^{X)), H : 
r*G' — ^ R, at any point X £ G defines a quadratic form on cotangent space T^G 
to G at X, and this quadratic form H\x '■= H\t^g is invariant with respect 
to the action of G on itself by left shifts. This fact is the motivation of the 
following definition. 

Definition 3.4. We say that a system with hamiltonian function p.40p de- 
scribes the free motion of a point on a Lie group G provided with left-invariant 
metric. 

Note that if the restriction H\x=e of H on T*G is positively defined, then H 
defines on G a left-invariant riemannian metric, and the system with hamiltonian 
H defines a geodesic flow on cotangent bundle T*G which corresponds to this 
metric. (H can be considered as the kinetic energy T of the point.) From 
a mathematical point of view, the positivity of the 2-form H\x=e, which is 
transferred on all G by left shifts, is not important (see TSI). 

Definition 3.5. Let the operator J in formula (I3.40p be diagonal, and let its 
element be such that 



with Ai > for i — 1, . . . ,n. In this particular case we also say that the 
system with hamiltonian function (j3.40l) describes the rotation of a rigid body 
in n-dimensional space, with generalized moments of inertia Ai, . . . , A„. In the 
following of this section we will always consider systems of this type. 

Lemma 3.13. If H = H{P^), then {H,P^} = 0. 

This result follows from relation (|3.10p . 



H = 




(3.40) 




(3.41) 
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Lemma 3.14. Let H = have the form {3.41^ , and suppose that Xi = Xj for 
some Then {H, P^} — 0. 

Proof. From (|3.8p and p.4ip we get for all z, j, I < i < j < n, 

n pL pL 

The right-hand side obviously vanishes when A.; — Xj. Note that the above 
formula is equivalent to the well-known Euler equations for the free rotation of 
a rigid body. □ 

Denoting A = (Ai, . . . , A„), let us group together the elements of the set A 
which are equal to each other. More precisely, we suppose that there exist u 
distinct positive real numbers /ii , . . . , /i„ , with 1 < w < n and 7^ Mfc for 
h ^ k, such that 

Ai = A2 = • • • = Apj = /ii , 

^pi+l = = ■ ■ ■ = Ap2 = ^2 , 

(3.42) 

+ 1 ~ -^p„-i+2 = • • • = Ap^ = , 

with pu = n. Let us consider the set q = q{X) = {qi, . . . ,qu), where qi = pi, 
q2 = P2 — Pi, ■ ■ ■ , qn = Pn — Pn-1- The set of all possible A can be divided into 
classes l{q) each characterized by a given set q = {qi, . . . ,qu), with qj e N for 
j = and X]J=i 1j ~ The order of the different generalized moments 

of inertia fih is not important, so one may always arrange them in such a way 
that 1 < 91 < 92 < • • • < (Zti- 

For a given A, consider the decomposition of R" which corresponds to this 

A: 

M" = Li e L2 © ■ • • ® , (3.43) 

where dimLj = qj, X G l{q), q = (gi, . . . , g„). Let us consider the subalgebra 
0(A) C = so{n) of all skew-symmetric matrices A such that the subspaces 
Lj are invariant with respect to the action of the operators defined by these 
matrices for all j = 1, . . . , g: 

9^ = {AeQ: A{L,) C Lj Vj = 1, . . . , u} . 

Let us consider a cartesian basis in R" corresponding to the decomposition 
p.43p . and the set of skew-symmetric matrices which represent the elements of 
Lie algebra 9 in this basis. A basis of g is then given by the set of N matrices 
of the form (12. 2p . with I < i < j < n. This basis is orthonormal with respect to 
the euclidean structure induced in q by the bilinear form {A, B) := Tr {AB)/2, 
for A, _B e g. Let denote the set of the pairs of indexes which correspond 
to equal generalized moments of inertia. Recalling equalities (I3.42p . we have 
^ (/!,...,/"), where 

/fc {(i,j) e N X N :pfe_i < i < j < pfe} for fc = l,...,w. (3.44) 

Then the set : (i, j) G of matrices of the form ()2.2p forms a basis 

of 0^. Let us consider the set P^^ := {Pi- : G I^} C P^ of impulses of 

the form p.7|) . The following proposition is an obvious consequence of lemmas 
m and EH 
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Proposition 3.15. If H = has the form then {H,B^} = 0, where 

Q\ := (pL^^pR). 

In other words, contains the elements of B which are in involution with 
the hamiltonian H^. We are now going to investigate whether it is possible to 
construct an integrable set using suitable functions of the set i?^ and possibly 
H"^. We will find that this depends on the properties of the set A or, more 
exactly, of the set of integers g(A). 

Let us consider a set of functions F = (Fi, . . . ,Fp) on the 2iV-dimensional 
symplectic manifold M"^^ . We suppose that rankF is the same almost ev- 
erywhere. Let k = k{F) be the maximal number of functionally independent 
functions, defined on all A'P^ , which are functions of F and are in involution 
with all functions of set F. More exactly, there exists a set Z = {Zi, . . . , Zk), 
Zj = Zj{F), such that these functions are functionally independent almost ev- 
erywhere on M and {Z, F} = on M, and k is the maximum integer for which 
such a set Z exists. 

Definition 3.6. We will say that the number k is the centrality of the set F, 
while r = r{F) :— 2N — rank F — k{F) is the defect of integrability of the set F. 

Lemma 3.16. For any set F (such that ranki^ is the same almost everywhere) 
we have k{F) < dimW^, where W :— Span di^H (Spandi^)^ at a typical point 
m € M . Moreover, r(F) > 0. 

Proof. For any function Z = Z{F), such that {Z,F} — at m, we must ob- 
viously have z e W. From this it easily follows that k{F) < dimM^. The 
inequality r{F) > follows instead from proposition 12. II □ 

Proposition 3.17. The centrality and defect of integrability of the set of func- 
tions B — {P^ , P^) are respectively 

KB) = [I] , (3.45) 

r{B)=Q. (3.46) 

Proof. Equality (|3.45l) follows from proposition 13.121 and remark 13.21 Then 
equality follows from ((X^ . □ 

Let : Q ^ denote the projector onto the subalgebra with respect to 
the euclidean structure in q introduced above. For any linear operator L : g — >• g, 
we denote by L\^\ its restriction to the subalgebra q^. For any m e M we can 
then consider the operators 

■P^ o adpi, : ^ 0^ , 
T'^oadpi.lgA :0-^^0^, 
adpL IgA : 0^^ ^ , 

where o denotes the composition of linear operators, and = P^{m). We 
define the three integers ai, and CT3, dependent on the point m, as the 
linear dimension of the kernels of the three above operators: 

:= dim Ker {V^ o adpi, ) , (3.47) 
diniKer (V^ o adpL\gx) , (3.48) 
CT^ := dimKer (adpi,|gA) . (3.49) 
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From these definitions and from p.25p it is immediate to see tliat > a > 
and CTi > CTj > 0-3. 

Proposition 3.18. The rank of the set = (P^^,P^) at any point m £ M 
is given by 

rankS^ = 2iV-cr^. (3.50) 

In addition we have 

dimW = cr + a^ - , (3.51) 

where 

Span dS^ n (Span dS^)^ . (3.52) 
Proof. Since = P^^ U P^, we liave 

(SpandP^)^ = (SpandP^)^ n (SpandP^^)^ 

= SpandP^ n (SpandP^^)^ , (3.53) 

where the last equality follows from p.l9p . Therefore (SpandP^)^ is the set of 
all those elements a e SpandP-^ such that Ii{a,dP^f.) = Oy{h,k) e I^. Any 
covector a € Span dP-^ can be expressed as a = J2i<j ^ij'^Pij ' where Aij are 
elements of a univocally determined skew-symmetric matrix A. Hence a = dP^ 
and 

n(a, dP^,) = {P^Ptk) = [P\AU , (3.54) 

where for the last equality use has been made of formula p. lip . The matrix 
elements on the right-hand side of the above formula, for (/i, k) G are just the 
components, with respect to the basis D'^ , of the projection P^([P^,v4]) of the 
element [P^jA] e onto the subalgebra q^. We see therefore that (Spandi?'^)^ 
is in one-to-one correspondence with the linear space of matrices A G so{n) such 
that P^([P^, A]) = 0. Hence 

dim(SpandP^)^ dimKer (P^ o adp^) = ct^ , (3.55) 

so that (|3.50p is obtained by applying lemma [3T1 to the set B^ . 
From and ([XS^ it follows that 

W = SpandP^ n (Span dP^)^ n (Span dP^^)^ . (3.56) 

Let us then consider an arbitrary element w €W . The condition w G SpandP'^ 
implies that there exist zi € Span dP^ and Z2 G Span dP^^ such that 

W = Zi + Z2 ■ 

Since P^^ C P^, from ((3T9)l it follows that 

H(zi,dP^^) = n(z2,dP^) = 0. (3.57) 

Hence the condition w G (SpandP''^)^ implies = H(w,dP^) = H(2i,dP^). 
Therefore, recalling p.l9p and (j3.20p . we have 

zi G Span dP^ n (Span dP^)^ = (SpandP)^ . 

Using proposition 13.71 we thus conclude that zi ~ J2i<j ^ij^^^ij' where C is a 
skew-symmetric matrix such that [P^,C] — 0, i.e., C G Keradpt. Finally, the 
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condition w S (SpandP-'^''')^, using again p.57p . implies — Il{w,dP'"'^) = 
U{z2,dP^^). Therefore 

Z2 e SpandP^^ n (SpandP^^)^ . 

Using formula (13.541) , we see that a covector a g Span dP^ belongs to Span dP^^C] 
(SpandP^^)-^ if and only if a = Y.^<] ^v^Pi^j^ with A e and V^{[P^,A]) = 
0. These two conditions on A can be simultaneously expressed as ^ € S^, where 
:=Ker (P^ o adp^ I^a) . 

We have thus shown that w — J2i<ji^ij + ^ij)dPij^ where C G 5* := 
Ker adp-L and A <^ . It follows that 

= |u; e T^M ■w = Y^ K^jdPtj , K e Span (S, S^)^ . (3.58) 

i<j 

Since covectors dPlj, 1 < i < j < n, are linearly independent (see proposi- 
tion [321), according to p.58p there exists a linear isomorphism between W and 
Span (5, 5^). Therefore 

dimVF = dimSpan(5,5-^) = dim 5* + dim S*^ - dim(5n 5-^) . (3.59) 

Recalling p.25p and p.48p . we have dimS* = cr and dimS*^ — Furthermore, 
it is easy to see that S (1 = Keradpt Ci = Keradpi|gA, so that recaUing 
([5^ we can write dim(5'n 5"^) = a^. Hence ([535)) is equivalent to ([XH]) . □ 

Corollary 3.19. At a typical point m £ AI, then the centrality k{B^) and 
the defect of integrahility r{B^) (see definition \3.6\) of the set satisfy the 
inequalities 

k{B^) <a + a^-a^, (3.60) 
r{B^) >(tI+(t^-g-g^. (3.61) 

Proof. Relation p.60p follows from lemma 13.161 and formula p.5ip . Then rela- 
tion ([33T|) follows from ([XSOl) and ([X60| . □ 

Note that formula (|3.52p implies dim < dim(Span rfP^)^. Therefore using 
p.55p and (|3.51l) we obtain the inequality > a + — cr^ ■ 

Let us consider again the set q = (gi, . . . , g^) introduced after formula (|3.42l) . 
If M = 1, i.e., q = (n), < Ai = A2 = • ■ ■ = A„, then = g, so that 
ai = a2 = cr^ = <J and B^ = B. We recall that, according to lemma 1331 we 
have a — [^] for a typical m e M. Therefore, in this particular case formulas 
([530)) and ([53g)) - ([53T)) agree with the resuhs rankP = 27V- [f ] , k{B) = [f ] , 
and r{B) ~ 0, see formula p. 321) and proposition 13. 171 

Lemma 3.20. Suppose that u > I, i.e., there exist at least two generalized 
moments of inertia which are different from each other. Then at a typical point 
m £ M we have 

1 / " \ 

(3.62) 



(3.63) 
(3.64) 
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On the right-hand side of formula i3.6S\) . the function d{q) is defined as the 
number of odd Qi, 1 < i < u or, equivalently: 



Proof. We shall outline the scheme of the proof for the case u — 2, i.e., q — 
((71,(72), with (71 + (72 = n. The generalization to the case of arbitrary u should 
be obvious. 

The decomposition (lOS)) of M" for m = 2 becomes M" = Li ® L2, with 
dimLi = (71, dimL2 = 92- According to this decomposition, the matrix 
P^{rn) G so(n) can be represented in a blockwise form as 



Ai B 

-B A2 



where Ai e so((7i), A2 G so{q2), whereas i? is a generic qi xq2 matrix. According 
to remark r3. 11 for a generic m d AI both matrices Ai and A2 will have pairwise 
different eigenvalues. By applying lemma l3^ to the two subspaces Li and L2, 
it is possible to find a basis in each of them such that both matrices Ai and 
A2 have the normal block-diagonal form p.2ip . In these coordinates, any other 
arbitrary matrix Z S so{n) can be represented as 

where Vi € so{qx), V2 € so{q2), whereas ?7 is a generic qi x (72 matrix. We have 
Z G ii and only ii U = 0. For the commutator of and Z we obtain 

[P\Z]^(% ;V (3.66) 



-D C2 



where 



Ci = [Ai,Vi] + UB- BU , (3.67) 
C2^[A2,V2]+UB- BU , (3.68) 
D = AiU -UA2-V1B + BV2. (3.69) 

We will have P^{[P^,Z]) = if and only if Ci = and C2 = 0. 

Let us now consider the number a2 defined by formula (|3.48p . The kernel of 
the operator V'^ o adpL \gx is made by the matrices Z of the form (j3.65p . with 
[/ = 0, such that Ci = and C2 = in ((3J6)) . Using formulas f^m\ - t^M^ . 
we find that Vi and V2 must satisfy the conditions [^i, Vi] = and [A2, V2] = 
respectively. According to formula p.22p we have that dimKerad^i = [qi/2] 
and dimKerad^2 = [(72/2], whence 



^'2 



'qi' 


+ 


"92" 


. 2 . 




. 2 . 



which corresponds to p.63p . 
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In a similar way, in order to evaluate the number tig defined by formula 
p.49p . we observe that the kernel of the operator adpi, |gA is made by the ma- 
trices Z of the form p.65p . with J7 = 0, such that [P^, Z] = 0. Using formulas 
p.66p - p.69p . we find that V\ and Vi must simultaneously satisfy the conditions 

[Ai,yi]=0, [A2,F2]=0, (3.70) 

and 

FiS-SV^2=0. (3.71) 

We recall that A\ and Ai are skew-symmetric matrices in the normal block- 
diagonal form (I3.21[) . with pairwise different eigenvalues. Hence, according to 
lemma conditions p.70p imply that also V\ and V2 must have the normal 
block-diagonal form (I3.2ip . with arbitrary eigenvalues. But then it is easy to 
verify that, for a generic i?, condition (I3.7ip necessarily implies that all eigen- 
values of V\ and Vi must be zero. From this we conclude that Z = 0, so that 
(I3.64P is proved. 

Finally, in order to evaluate the number a\ defined by formula p.47p . we 
note that the kernel of the operator o adpi, is made by the matrices Z of the 
form (j3.65p . such that Ci = and C2 = in p.66p . Introducing the matrices 
Ti BU - UB and T2 BU - UB, for such Z formulas (|3171)-(I31S1) can be 
rewritten as 

[Ai,Vi]^Ti, [A2,V2]=T2. (3.72) 
Hence U must be such that 

Ti e adAi (01) , T2 e adA, (02) , (3.73) 

where gi — so{qi) and 02 — so{q2). According to formula p.24p . this is equiva- 
lent to 

{Tl)l2 = (Ti)34 = • • • = (ri)2si-l,2si = , 
(72)12 = (72)34 = • • • = (72)2^2-1, 2S2 = , 

where si = [qi/2], S2 — [92/2]. It is easy to check that, for a typical matrix 
B, the linear system of equations (I3.74p can be solved with respect to si -|- S2 
appropriately chosen elements of the matrix U. Therefore matrices U satisfying 
conditions p.73p form a linear space having dimension qiq2 — si — 82- After 
choosing U in this space, in order to obtain an element Z G Ker (V^ o adpi,) 
one has to take Vi and V2 satisfying equations p. 721) . These are equivalent to 

Viead^liT,), V2ead^l{T2). 

According to coroUarv 13.51 we have 

dimad^^(Ti) = dimKerad^ij — si , 

dimad^^(T2) — dimKerad^2 — S2 ■ 

More exactly, equations (I3.72p determine all elements of Vi and V2 except (Vi)i2, 
(^"1)34, • ■ • , (Vi)2si-i,2si, and (^2)12, (^2)34, • ■ • , (V2)2s2-i,2s2- These are S1+S2 
elements which can be arbitrarily chosen, in addition to qiq2 — si — S2 elements 
of the matrix U previously considered. We thus conclude that 

O"! = gi'72 , 

which corresponds to p.62p . □ 
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Let denote the standard set of Casimir functions for the adjoint action 
of the group SO{qj) of orthogonal transformations on euchdean subspace Lj C 
R", j — l,...,u, see definition 13.21 This set contains [(?j/2] functions on the 
algebra so{qj) which corresponds to the subspace Lj C R". We can consider 
these functions as functions on the algebra g — so{n). Let us consider the set 

— {C-^, . . . , C") of functions on g, obtained by collection of sets C^ . Clearly 

contains 

u 



^ (3.75) 



elements. Let C^^ = . . . , C^") denote the set C^^ := C^oP^ of functions 

on r*G, obtained by making the composition of the functions of set with 
the map . We have already introduced the set C = C o P^ = C o P^ (see 
proposition l3.1ip . Clearly C^'^ = C if u = 1. We will denote with Z'^ the set of 
functions := {C, C^^), if u > 1, or Z^ C, if u = 1. The set Z^ contains 



elements, where 




z^^v;';\ , (3.76) 



Proposition 3.21. For any X — (Ai, . . . , A„), each function of the set Z^ is in 

involution with each function of the set : 

{Z^,B^} = Q. (3.77) 

Furthermore, almost everywhere in M = T*G, where G — SO{n), the set 
Z^ is functionally independent, i.e., i&nkZ^ ~ z^. 

Proof The equality {C,B^} = follows from ((336| . From the definition 
of Casimir function, and the fact that is a Poisson map, it follows that 
^(jLi pLzy = oVi = where P^^ := {P^^^ : {h,k) € P}, see formula 

(EUl). Relation dXH]) implies {P^\P^3} = OVi ^ j, whence {C^\P^^} = 0. 
Furthermore, (IXTU|) implies {C^^,P"} = 0. Hence {C^^,B^} ^ 0, and equal- 
ity (j3.77p is proved. 

According to proposition l3.12l applied to the group SO(n) and its subgroups 
SO{qj), j — I, . . . ,u, both sets C and C^'^ are almost everywhere functionally 
independent, so that 



rankC^^ = s^. (3.78) 



rank C = 

It follows that 

rank Z^ = dim Span {dC, dC^^) 

= I +s^-dim(SpandC'nSpandC^^). (3.79) 

Let us then consider the set 

G = |w e r;„Af ■■w = Y, ^^j'^Ptj ,K£Ker (adpi, |j,a) } . (3.80) 

Using ((X^ and we see that 

G = Span dCn Span dP^^. (3.81) 
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Since SpandC^'*' C Span dP-^-** , we obtain 

GD Span dCn Span dC^^ 



(3.82) 



On the other hand, any w G G can be expressed as w = X^iLi "^ith Wi e 
SpandP-^*. For any z G SpandP-^* we have — Il{w,z) — Il{wi,z). Hence 
Wi e SpandP-^* n (Span c?P^')^ = SpandC-^*, where the last equaUty follows 
from the application of (j3.39p to the subspace L^. It follows that w G Span dC^^. 
Hence, using (|3.8ip and p.82p . we conclude that 



G = Span dCn Span dC 



LA 



SO that 



dim (SpandC n SpandC ) = dimG = dim(Keradpi,|jjA) = 



(3.83) 



Recalling formula (j3.64p . and the equalities s''^ — — [n/2] for u = 1, we 



then obtain from p.79p that rankZ"^ 
functionally independent. 



for any u, so that the set Z is 

□ 



Lemma 3.22. Almost everywhere in T*G we have 

SpandZ^ = SpandP^ n (SpandP^)-^ 
Moreover, if u > 1 we have 



(3.84) 

(3.85) 
(3.86) 
(3.87) 



where k{B^) and r{B^) are respectively the centrality and the defect of integra- 
bility (see definition \3.b]) of the set . 




V 


\d{q) + 11 
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Proof. Since the elements of are functions of B^, from equality (13.771) it 
follows that SpandZ^ C W -.^ Spa^ndB^ D (SpandS^)^. From (I53T|) . taking 
into account p.27p . (|3.63l) and (|3.64p . we obtain that dimW^ — z^. Moreover, 
we know from proposition 13.211 that the set Z^ is functionally independent, so 
that dim(Span dZ'^) = z^. Hence we conclude that SpandZ^ = W, so that 
is proved. 

Equality (p:S5)) follows from (P^Ul) and (IX^ . 

According to definition 13.61 the centrality of B^ is at least as great as the 
number of elements of Z"^, i.e., we have k{B^) > z^. On the other hand, using 
(1323), (I5:M)) and (13211), we obtain from ^^W^ that fc(P^) < z^. Therefore 
(I3.86P is proved. 

Finally, p.87p follows from definition 13.61 and from equalities p.85l) - p.86p . 

□ 
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Note that one can easily prove, by induction on u, that r{B^) given by 
formula p.87p is always an even integer. 

The complete integrability of the system describing the free rotation of an 
n-dimensional rigid body can be proved by introducing the so-called Manakov's 
integrals [20]. It is not difficult to check that the function = (l/2A:)Tr (P^ + 
J^pY ■ T*G — >■ R, where J is the diagonal n x n matrix with Ai, A2, . . . , A„ as 
diagonal elements, is in involution with i?'*' for any value of the parameter 
p. Hence the coefficients of the polynomial 'P{p) in the variable p are also in 
involution with H^, i.e. we have {cij,H'^} = 0, where 

k 

V{p)^^TT{P^ + J^pf=Y,cu,p>. 

These coefficients are not all functionally independent. One immediately sees 
that Cfcfc is just a constant, and that c^j = whenever k — j is odd. Furthermore 
it can be proved that, if one is only interested to functionally independent 
elements, then one need only consider coefficients Ckj with k — 2,S, . . . ,n. It is 
easy to see that 'P{p) is in involution with all functions of the set B^, so that 
{cij,B^} — 0. Moreover, one can prove that all coefficients are mutually in 
involution, {cij,Ciiji} = 0. This result provides in particular another proof of 
the fact that all these coefficients are integrals of the system with hamiltonian 
H'^, for it can be shown that i?^ can be expressed as a linear combination of 
the functions Ck,k-2 for fc = 2, . . . , n. 

It has been proved in general that the system with hamiltonian is in- 
tegrable [211 US HZl ■ In the general case in which all generalized moments of 
inertia are pairwise different, A^ Xj for i 7^ j, an integrable set of functions is 
given by {M]P'^), where M = {ck,k-2i, k = 2,...,n, i = 1, 2, . . . , [fc/2]) is the 
complete set of X]fc=2['''/-^] ~ {^/'^){^{^^ + ["-/S]) functionally independent 
coefficients c^, and P'^ is a set of n{n — l)/2 — [n/2] elements of P^, such 
that {C, P'^) is a functionally independent set. Hence this system is integrable 
with (l/2)(n(n — l)/2 + [n/2]) central functions. The elements of M such 
that J > are called Manakov's integrals. The remaining [n/2] elements of 
M, i.e. C2kfi = (l/2fc)Tr(P^)'^' for k = 1, 2, . . . , [n/2], are independent of the 
moments of inertia A, and form a set of Casimir functions equivalent to the set 
C introduced in section 13.11 Hence, an equivalent integral set of functions for 
the free n-dimensional rigid body with pairwise different moments of inertia is 
(C, M; P'^), where M is the set of (l/2)(n(n-l)/2- [n/2]) Manakov's integrals. 

When the moments of inertia are not all pairwise different, the set M is 
no longer functionally independent. However the integrability of the system is 
preserved, which means that one can construct an integrable set of functions 
whose central subset is made of the elements of and of a suitable subset 
of M . According to proposition 12.11 such a subset of M must contain just r/2 
elements, where r — r{B^) is the defect of integrability of the set B^ and is given 
by formula (|3.87l) . Therefore, the central subset will contain k{B^) +r{B^)/2 
elements. This result is expressed by the following proposition. 

Proposition 3.23. The system with hamiltonian H = given by formula 
US. is integrable with 

Hq)^\(^n' + 2n-j2lf^-^^ 



d{q) + 1 



(3.88) 
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central integrals. 

Manakov's integrals can be explicitly represented in the following form: 

Ck,k-2l = ^ "-Vk-T Pki^P^ts ■ ■ ■ PL-ii'ii^Ln ' (3-89) 

11,12, ■■■,121 

with < ^ < fc/2, where 

<r2r- E \t^xt^---\t:^5,,+,,+...+,,,.k-2i. (3.90) 

6i>0,b2>0,..^,fc2!>0 

We see that Ck.k-21 is a homogeneous polynomial of degree 21 in the left-invariant 
momenta, while its coefficients a^^Jj'll2z' homogeneous polynomials of degree 
2(/c — 21) in the generalized moments of inertia, completely symmetrical with 
respect to permutations of the indexes ii, ... ,121- 

In Table [3] we give the number k{q) of central integrals resulting from the 
above proposition for free n-dimensional rigid bodies with n < 6. We also give 
the quantities k{B^) and r{B^) resulting from lemma [3.221 



n 


q 


k{B^) 


r(B^) 


k{q) 


3 


(3) 


1 





1 


3 


(1,2) 


2 





2 


3 


(1,1,1) 


1 


2 


2 


4 


(4) 


2 





2 


4 


(1,3) 


3 





3 


4 


(2,2) 


4 





4 


4 


(1,1,2) 


3 


2 


4 


4 


(1,1,1,1) 


2 


6 


5 


5 


(5) 


2 





2 


5 


(1,4) 


4 





4 


5 


(2,3) 


4 


2 


5 


5 


(1,1,3) 


3 


4 


5 


5 


(1,2,2) 


4 


4 


6 


5 


(1,1,1,2) 


3 


6 


6 


5 


(1,1,1,1,1) 


2 


8 


6 


6 


(6) 


3 





3 


6 


(1,5) 


5 





5 


6 


(2,4) 


6 


2 


7 


6 


(3,3) 


6 


2 


7 


6 


(1,1,4) 


5 


4 


7 


6 


(1,2,3) 


5 


6 


8 


6 


(1,1,1,3) 


4 


8 


8 


6 


(1,1,2,2) 


5 


8 


9 


6 


(1,1,1,1,2) 


4 


10 


9 


6 


(1,1,1,1,1,1) 


3 


12 


9 



Table 3: Number of central integrals for free rigid bodies. 
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3.3 Free rotation of a quantum rigid body 

In order to quantize a free rigid body we have to consider the quantum impulses 
i^' ii' '^tiich constructed according to formula (j2.16l) in correspondence 
with vector fields V^^, respectively, 1 < i < j < n. However, if we want 
to apply to this system the concept of integral quantum system introduced in 
[T], we are apparently faced by the problem that here the configuration space 
K — G — SO{n) is not a domain of the linear space M.^ . This problem is 
solved by the consideration of local coordinates on G. Note that the main 
part MT of a linear differential operator (see definition in 1 ) is not defined 
intrinsically. However the symbol (MJ^)*^™*^ can be considered as intrinsically 
defined, according to the following proposition. 

Proposition 3.24. The symbol S := {MTf^^ of the main part of a linear 
operator of class O, expressed via local coordinates x on configuration space 
K , has the form of a homogeneous polynomial of p. This polynomial S — 
S{x,p) behaves under a change of local coordinates x on K as a function on 
the cotangent bundle T*K to the manifold K. It follows from this fact that 
the definition of quasi-independence of a set of operators does not depend on 
the choice of local coordinates on configuration space K . The same is true for 
the definition of quasi-integrability of either a set of operators or an individual 
operator. 

The proof of this proposition is obvious. The first part of the proposition, 
about the representation of S* as a function on T*K, is actually the reformulation 
of well-known facts. 

Let us consider the quantum system with hamiltonian operator 



H = 



- , 2 



^EF^ (3-91) 



Xi + Xj 



on G°°{SO{n)). We consider this system as the system describing the free 
rotation of a quantum n-dimensional rigid body. 

Proposition 3.25. For n < 6 this quantum system is quasi-integrable for any 
X, with the same number k of central operators as the number k{u) of central 
integrals of the corresponding classical system, see proposition \ 3. 23\ and Table 

m 

Moreover, if q = {n), this quantum system is quasi-integrable for any n with 
[n/2] central operators. If q = (l,7i— 1), this quantum system is quasi-integrable 
for any n with n — 1 central operators. 

Proof Let and P^ denote the sets of operators P^ = {P^ , 1 < « < 
j < n) and P^ — {PijA < * < 3 < "-)• Let us consider the set B :— 
{P^,P^) containing 2N operators. Let denote the set of operators which 
are obtained by symmetrization with respect to P^ from the functions of set 
G^ = C(P^), i.e., := (G^yy'^. Analogously we define C^, P^^, B^, 
C^^, Z^, i.e., := (C^f^"' etc. We first show that the commutation re- 
lations [13^, H^] = [B^,Z^] = follow from the analo gous classical relations 
{B^, H^} = {B^, Z^} = and from the propositions about quantization of [5] 
and section [32] of the present paper. The commutators of the operators of set 
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B have the same form as the Poisson brackets of the corresponding classical 
functions, since all these functions are linearly dependent on canonical impulses 
p, see formulas p.7l) - p.l0p and proposition 12.61 

[^^ , PL] = -S.^Pji - ^jkPth + 5^kP^h + hhPk , (3-92) 
[Af , P« ] = S^^Pf^ + S.kPi - S^kPjl - S.hPi , (3.93) 
[A^,P«]=0. (3.94) 

The relation [H^,B^] = then follows from {H^,B''^} = using proposition 
4.2, case b, of [5]. Similarly, since 

{C^(F^),B^} = 0, = (C^)«>"° , 

the equality [C^,B'^] = is obtained by applying corollary 4.4 of [2]. Analo- 
gously, we obtain the equality [C^^,B'^] — from the corresponding classical 
relation {C^^, B^} = 0. We have thus proved that [B^, Z^] = 0. 

liq = (n) 01 q — (1, n—1), according to lemma [?.22l the defect of integrability 
of the set B^ in the classical case is r{B^) — 0. Moreover, we have k{B^) — [n/2] 
if q = (n), and k{B^) — n ~ 1 ii q = {l,n — 1). This implies that in these 
two cases, for any n, there exists a classical integrable set of functions of the 
form F = {Z^; B'), where B' C B^, and the central subset Z-^ contains k{B^) 
elements. Let us then consider the corresponding set of operators F — {Z^; B'). 
From what we have seen above, it follows that [Z^, F] = 0. Moreover, since all 
functions of F are homogeneous with respect to p, these functions coincide with 
their main parts with respect to p, i.e., M{F) — F. It is also easy to see that 
the elements of F are the symbols of the main parts with respect to p of the 
elements of the corresponding set of operators F, i.e., F = {MPy^^ . Hence, 
the quasi-independence of the sets of operators F follows immediately from the 
functional independence of set of functions F. One thus concludes that F is 
an integrable set of operators with k{B^) central elements. One can also easily 
show that in these two cases, H'^ is a linear combination of the elements of Z^. 
The integrability of the system describing the free quantum rigid-body is thus 
proved for any n in the two cases q = (n) and g = (1, n — 1). 

In the remaining cases, the classical integrable sets of functions generally 
include also one or more Manakov's integrals among their central elements. 
We define Manakov's operators Ck,k-2i as the symmetrization of the classical 
functions p.89p with respect to the left-invariant momenta: 

£k.k-2l = li E <!^-2rSyr^2liPn^.,PL,■■■,Pt-^^..,PLJ^ (3-95) 

11,12, ■■■,i2l 

with < I < k/2. Note that, also in the quantum case, the hamiltonian operator 
H'^ can be expressed as a linear combination of the operators Ck,k-2 for k — 
2, . . . , n. By applying again the results of ^2, we easily see that [ck,k-2iT B'^] = 0. 
However, no general theorem ensures that Manakov's operators commute with 
or among themselves. We will here limit ourselves to studying commutators 
between Manakov's operators of degree lower that 6 in the momenta. This will 
be sufficient to establish the quasi-integrability of the free quantum rigid body 
in spatial dimensions n < 6. 
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The commutator between two Manakov's operators, when one of them is 
of second degree, can be evahiated by making use of proposition 2.3 of [5], 
and of the algebra p. 921) of left-invariant momenta. In this way, with some 
computation we find for any I, h: 

[ci,i-2,Ch,h-2\ = 0, (3.96) 

[ci,i-2,CH,k-4] = IY. Kh^y-^^iPtv P^i.Pti) , (3.97) 

where 

lijk _ ij ( r) iijk o.,iikk \ ' iikp \ 

+ E <r2(<t4-<^4)- (3.98) 

Note that the operator SyiR^{P^ , Pj^j^, P^^) is completely antisymmetrical with 
respect to permutations of indexes i, j, fc, so that formula p.97p can be rewritten 
as 

[q,/_2,c,,„_4]= E fcff Sym3(i^^,i^iP^,), 

i<j<k 

where b'^ij^^ denotes the complete antisymmetrization of coefficient b^^ . 

In formula p.98p the coefficients a;'';_2 and al^f^_/^ have to be replaced by 
their explicit expressions given by (I3.90p . We have 

1-2 x2(/-l) _ \2(/-l) 

4-2 = E A?('--)Af = . (3.99) 

Moreover, for h — 5 we have 

i]kp _ ijkp _ x2 I x2 I 1,2 I ^2 
^hJi-4 — "5,1 — -^i + ^k + \ ^ 

SO that from p.98p we easily obtain b^i^^^ = 0, which implies 

[Q,i_2,C5,l] =0 (3.100) 

and consequently [H'^,c^^i\ = 0. From these results it follows that the free 
quantum rigid body is a quasi-integrable system for spatial dimensions n < 5. 
The integrable set of functions of the classical system can in fact be quantized 
by replacing Manakov's integrals with the corresponding operators. Since these 
operators are homogeneous in the momenta, the symbols of their main parts 
with respect to p coincide with the corresponding classical functions. Hence the 
quasi-independence of the sets of operators is a consequence of the functional 
independence of the classical sets of functions. 
For h — 6 the coefficient aj^^f^_^ becomes 

^yfep — )i \i _|_ ^4 ^2y2 ^2^2 _|_ ;\2^2 _|_ ^2^2 _|_ \^.\^ A^Ap , 
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and from p.98p we obtain 



_ 5 
~ 6 



A^(A^ - + Apa0_2 + >^i^ A^)«u-2 (3.101) 

Xf-'' (A^ - AD + Af (A^ - A?) + Af (A? - A^) 



Since for a generic set A of generalized moments of inertia the above expression is 
different from 0, we have that in general [q_/_2, C6,2] 7^ 0. Note also that, putting 
I = 3/2 in (I3.99p . we get a^^^ ~ l/{Xi + Xj), so that we can formally write 
= — C3/2,-i/2- From the above formulas we thus directly obtain 



i<j<k 



ki) 7 



with 



5U* 



5 r 
6 



A.(A2-A2) + A,(A2_An + Afe(A2-A?) 



Hence C6,2 is not in general an integral operator of the quantum system with 
hamiltonian H^. 

By using again proposition 2.3 of 2 and commutation relations p. 921) . one 
finds however that, for arbitrary symmetrical coefficients a^^ = a- 
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4 s: 



i<j<k 



with 



a''')«U-2 + (a 



a-' )a-' 



Ll-2 



1,1-2 



(3.102) 



One sees immediately that the above expression becomes identical with p.lOip 
if one chooses a'^ — (5/6)Af A^. Hence, if we define the modified Manakov's 
operator 



6,2 



^6,2+^y:A?A|(^^)^ 



12 ^ 



then we get 



[Cl.l- 



-2, 6*6,2] — , 



66,2] = . 



For n = 6, let us then consider the set of operators F which is obtained from 
the classical integrable set of functions F given by proposition l3.23l by replacing 
Manakov's integrals c/,i_2 = 3, . . . , 6), 05,1, and C6,2, with the operators c;,/_2, 
£5,1, and C6,2 respectively. In order to prove that this set of operators satisfies 
the required commutation relations, we still have to show that 



[55,1,66,2] -0. 
With the techniques already employed one obtains that 



(3.103) 



i|Ea?a?[ 



T^> ■A?A^fc5,i,(f^)2] 



\4\2 



6^' 

hlr, 



AfA^ -Sym3(P^;, P;,„, P^^j^ 



E ^y™5 i^tj 1 Pjh ' Phi ' Plm ' Pmi) 



(3.104) 
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On the other hand, the commutator [£5^1, £6,2] between two fourth-order sym- 
metrized polynomials cannot be worked out with the tools provided in [2]. By 
means of a straightforward and quite heavy calculation, we have verified that 
this commutator is indeed just the opposite of the expression (I3.104p . so that 
equality (|3.103p actually holds. Of course, the symbol of the main part of (^6,2 
coincides with the classical function C6.2, so that the quasi-independence of the 
set F again follows from the functional independence of the classical set F. We 
conclude that -F is a quasi-integrable set of operators. The proposition is thus 
completely proved. □ 

We have seen that, for n > 6, the correct quantization of Manakov's integral 
C6.2 does not coincide with the symmetrization of the classical function with re- 
spect to the left-invariant momenta. We have nevertheless provided a recipe to 
achieve the quasi-integrability of the free quantum rigid body for n = 6 for ar- 
bitrary moments of inertia A. We can conjecture that analogous procedures can 
lead to the quasi-integrability of this quantum system also for n > 6. However, 
we are unable at the moment to prove that this conjecture is true. 
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